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ABSTRACT 
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Control/Life Support System Computer 
Program 

by 

Ernest E, Smith, Jr. 


Contract NAS 9-14913 


February 1978 


This users guide describes the computer programs developed to aid 
in the analysis of RSECS (Representative Shuttle Environmental 
Control System) related material. These programs have been pre- 
pared to provide pretest predictions, post-test analysis, and real- 
time problem analysis for RSECS related test planning and evalua- 
tion. The programs have been designed for use on a Wang 2200 series 
computer system. Hamilton Standard has provided these programs to 
the NASA on magnetic tape cassette cartridges and on a disk device 
that is part of Crew Systems Division's Wang 2200 series computer 
system. 
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the requirements of Contract NAS 9-14913, Shuttle/RSECS Computer 
Program. This report covers work accomplished during calendar years 
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INTRODUCTION 


To fulfill the requirements of Contract NAS 9-14913, for calendar 
year 1977, Hamilton Standard has developed the computer programs 
listed below. These programs were written to support RSECS related 
analysis. 

- "SUITSK2" - Calculates the steady state parameters of the EMU 
(Extravehicular Mobility Unit) for the IVA (Intravehicular 
Activity) mode when it interfaces with the Shuttle environmental 
control system. The program also calculates the EMU parameters 
in the EVA (Extravehicular Activity) mode. Because of its size, 
the program does not provide for a printout of the results. 

This function is provided by program "SUITOUT". 

- "SUITOUT" - Provides the printout of the parameters calculated 
in "SUITSK2" program. The printout is provided in a tabularized 
format and an optional flowchart format. 

"LIOHPPC" - Plots the partial pressure of CO 2 at the upper and 
lower LiOH canister outlets and the ARS inlet versus time. The 
partial pressures of CO 2 were provided by the RSECS LiOH test 
data. The program stores the input data on cassette tapes and 
provides an optional printout of the data in a tabularized format. 

- "LIOHCRT" - Plots the instantaneous removal rate of CO 2 from the 
LiOH canisters. The program produces a separate plot for the 
upper and the lower canisters using the RSECS LiOH test data pre- 
viously stored on cassette tapes by the "LIOHPPC" program. 

- "LIOHMAC" - Plots the accumulation of CO 2 in the LiOH canisters 
versus time. The program produces a separate plot for the upper 
and lower LiOH canisters using the RSECS LiOH test data previously 
stored on cassette tapes by the "LIOHPPC" program. 

- "BLOCKF” - Numerically approximates the radiation view factor be- 
tween a rectangle and an N-sided p&lygon. The program operates 
on two surfaces at a time and will not include the effects of a 
shadowing surface. 

"SCRIPTF" - Calculates radiation exchange factors from a user 
provided view factor matrix, area vector, and emitance vector. 


SUITSK2 STEADY STATE COMPUTER PROGRAM 


File Name "SUITSK2" 

Abstract "SUITSK2" calculates the steady state temperatures, 

pressures, flow rates, and heat loads around both the 
air and water loops of the EMU (Extravehicular Mo- 
bility Unit) for both the IVA (Intravehicular Activity) 
and the EVA (Extravehicular Activity) modes. Steady 
state solutions can be achieved for a range of meta- 
bolic rates and heat leaks which simulate anticipated 
operational requirements. The program is designed for 
use with the Wang 2200 series computer system. 

Program Description 

The program was written in Fortran by Hamilton Standard to assist in 
EMU performance analysis. The program has been converted into Basic 
and split into two separate programs for use on the Wang 2200 system. 
The first of the two Wang programs ("SUITSK2) requests input data 
and performs calculations, the second ("SUITOUT") prints out the 
results. The original Hamilton Standard program has been through 
several revisions. The version of the program which was translated 
into Basic was " SUITSK2" which was the latest available version. 

The program typically works in the "no sweating" mode. Skin temper- 
ature is a function of total metabolic rate. If the system cannot 
meet the heat removal requirements, a sweating mode is initiated. 

The maximum sweating rate is determined by the maximum relative hu- 
midity leaving the suit (an input value). If the heat removal rate 
can still not be met with sweating, the program reduces the total 
metabolic heat load until the system can handle the heat load. The 
difference in total metabolic heat loads must be handled by the man 
as stored heat. 

It should be noted that only one skin temperature can occur at a 
given metabolic rate. In actuality, a range of skin temperatures are 
possible at each metabolic rate. The program is presently using the 
maximum skin temperature (sweating threshold) vs. metabolic rate. 

The program attempts to fix the LCG heat load from the man as a func- 
tion of metabolic work rate approximating available NASA test data. 
The balance of the heat load is then split between the suit sensible 
and suit latent heat loads. At higher suit pressures, the sensible 
heat load is so high that QLCG must deviate from the curve in order 
to maintain a heat balance. 

The water loop heat exchanger, the suit, and sublimator when used as 
a heat exchanger, are all handled as counterflow heat exchangers 
(Q = UA x DTLM) . The sublimator is actually a cross flow heat ex- 
changer but due to the high mass flow ratio (Wcp)H 20 /W C p)AIR) , it 
can be considered to be counterflow. 

The program adjusts sublimator water loop conductances for sublimator 
water loop flow rate and average temperature. 


The program converges on two variables and must go a minimum of four 
iterations. The convergence variables are: 

1) Suit inlet gas temperature 

2) Skin temperature 

If the above variables are within tolerance for two successive itera- 
tions, the program is considered converged. 

Two special functions, DEFFN1 and DEFFN2, are included in the program. 
DEFFN1 calculates a saturation pressure in psia based upon a saturation 
temperature in degrees F. The function was derived from a sixth order 
regression of the steam tables. DEFFN2 calculates a saturation temper- 
ature in degrees F based upon the log of 10 times the saturation pres- 
sure in psia. The function is based upon a subroutine "KANDK" used by 
the original Fortran version of the program. 

The input data is identified in Table I and is requested by the CRT at 
the beginning of the program execution. Table I also identifies typi- 
cal input data for a variety of runs. Before actual calculations begin, 
the input data is printed out so that it can be verified. An example 
of the input data printout is shown in Figure 1. 

Next, the CRT asks whether an optional diagnostic printout is desired. 
The diagnostic printout provides an incremental printout of the varia- 
bles as they are calculated. Otherwise, interum variable calculations 
are not printed out. The diagnostic printout is usually not required, 
however, when difficulties occur during a run, a rerun with a diagnos- 
tic printout will aid in identifying the problem. 

The program iterates until either an error occurs, the maximum number 
of iterations is reached (50 iterations), or the convergence criteria 
is met. The first two conditions will result in an error message on 
the CRT, while the latter will result in a loops converged message on 
the CRT. At this time, the CRT will instruct the user to load the se- 
cond program ("SUITOUT") A listing of "SUITSK2" is provided in 
Table II. Table III provides a cross reference between the original 
Fortran variable names and the Basic variable names in the event that 
further program changes are required. 


Figure 1 


SUITSK2 PROGRAM EXAMPLE OUTPUT 


SUITSK2 PROGRAM RUN 


PROGRAM INPUT DATA: 

G.METT, GENV, QFAN, QPUMP= 1000 0 11 145 

PTOT, PC02, CFM, WH20, UAST- 15.2 9. 67000000E-03 6.6 185 60 

TWI , WCPW, CLIOH, CL2~ 45 600 6 3.2 

CL3, CL4, RH= 15 50 .9 

TFI.TSIJBO, XKEYJJAX, WLIOH, AP,UA.W= 115 50 1 864.6 21 .2361 601 

RFIN. HAG AS, GDP, IJAHEL.= 471.6 16.28 5. 42000000E-03 .4 

TSKINT, QDCMEN, QDCMEL, QEVCS, QPLSEN, INC= 95.25 0 27 41 0 1 




TABLE 1 SUITSK2 INPUT DATA DEFINITION 


NAME 

DESCRIPTION 


TYPICAL VALUES 

UNITS 

QMETT 

TOTAL METABOLIC LOAD 


UI* 

BTU/HR 

QENV 

ENVIRONMENTAL HEAT LOAD TO SUIT 


UI 

BTU/HR 

QFAN 

FAN POWER TO AIRSTREAM 


UI 

BTU/HR 

QPUMP 

PUMP POWER HEAT LOAD 


UI 

BTU/HR 

ptot” 

TOTAL PRESSURE - GAS 


UI 

PSIA 

PC02 

-C02 PARTIAL PRESSURE INTO SUIT 


.00967 

PSIA 

CFM 

FAN VOLUME FLOW RATE 


UI 

FT 3 /INCH 

WH20 

PUMP FLOW RATE (H20) 


UI 

LBS/HR 

UAST 

UA OF LCG (H20) 


UI 

BTU/HR-F 

TWI 

COOLANT H20 INLET TEMPERATURE TO SUBLIMATOR 


45 

F 

WCPW 

COOLANT H20 MASS FLOW RATE IN VEHICLE HX 


600 

BTU/HR-F 

CLIOH 

CONDUCTANCE BETWEEN LIOH AND CANISTER 


6 

BTU/HR-F 

CL2 

CONDUCTANCE BETWEEN GAS AND CANISTER 


3.2 

BTU/HR-F 

CL3 

CONDUCTANCE BETWEEN CANISTER AND WATER LOOP 


15 

BTU/HR-F 

CL4 

CONDUCTANCE BETWEEN LIOH AND GAS STREAM 


50 

BTU/HR-F 

RH 

MAXIMUM RELATIVE HUMIDITY LEAVING SUIT 


.9 

— 

TFI 

GUESS OF FAN INLET TEMPERATURE 

o © 

115 

F 

TSUBO 

GUESS OF GAS SUIT INLET TEMPERATURE 

50 

50 

F 

XKEY 

TYPE OF SYSTEM - 1=IVA, 2=EVA 


UI 

— 

UAX 

UA OF VEHICLE HX 


864.6 

BTU/HR-F 

WLIOH 

H20 FLOW TO LIOH CANISTER FOR COOLING 


UI 

LBS /HR 

AP 

POROUS PLATE AREA 


.2361 

FT 2 

UAW 

UA OF H20 LOOP IN SUBLIMATOR 


UI 

BTU/HR-F 

RFIN 

SUBLIMATOR FIN RESISTANCE 


471.6 

BTU/HR-F 

HAGAS 

UA OF GAS LOOP IN SUBLIMATOR 


UI 

BTU/HR-F 

CDP 

PRESSURE LOSS BETWEEN SUIT INLET TO FAN INLET 

K to 

.00542 

PSI/ (LBS/HR) 1 - 3 

UAHEL 

HELMET CONDUCTANCE TO INLET GAS STREAM 


.4 

BTU/HR-F 

TSKINT 

INITIAL SKIN TEMP AT REST 


UI 

F 

QDCMEN 

ENVIRONMENTAL HEAT LOAD TO DCM 


UI 

BTU/HR 

QDCMEL 

ELECTRICAL HEAT LOAD TO DCM 


UI 

BTU/HR 

QEVCS 

ELECTRICAL HEAT LOAD TO EVCS 


UI 

BTU/HR 

QPLSEN 

ENVIRONMENTAL HEAT LOAD TO PLSS 


UI 

BTU/HR 

INC 

INC-FACTOR TO ALTER SUBLIMATOR UA (TYPICALLY EQUAL TO 1) 

UI 

— 


*UI-user input. The user must define required values. 


TABLE II SUITSK2 PROGRAM LISTING 


10 REM SUITSK2 PROGRAM 

20 COM 01(53), Q2(18),Q3(19),G4(8),Q5(8>, 06(14) 

30 COM U1 ( 4 ) , A1 , A2, B , C , Cl , C2, C3, C4, C5, D, D1 , D7, El , E2, G1 , HI , K 
40 COM PI (2 ) , P2( 16 ) , P3 ( 7 ) , P4 (6 ) , 06, OS, Rl, R5, R6, S, Tli R4 
50 COM W1 ( 19) , W2( 1 ) , W3 ( 1 ) , U1 , U2, U3, U4, U5, U6, U7, US, U9, V2, V5, W4, W5 
60 COM T1 ( 24 ) ,T2(7) ,T3( 12) ,T4( 29 ) , T5(2) , T6( 2) 

70 COM VI ( 17 ) , V2( 8 ) , V3( 16 ) , V4 ( 7 ) , I 
80 DIM Y(4,4),Z(4, 1),X(4, 1), 1(4,4) 

90 DEPFN1 ( X ) -2. 78076237E-02+5. 13243277E~04*X+3. 98156450E-05*Xt2- 
4. 66483844E-08*Xt3+4. 91065757E-09*XT4- 1 . 779541 14E-12*Xt5+4. 63530 
688E” 14*Xt6 


100 DEFFN2 ( X ) = 35. 1 5789+24. 592588*X+2. 1 182069*Xf 2- . 3414474*X*3+. 1 
5741642-K-Xt4-„0313295S5*-Xf5+,.0038658282*Xt6-.00024901784*X1’7+.QOO 
Q068401559*Xt8 

110 PRINT "TOTAL METABOLIC LOAD": INPUT 01(8) 

"ENVIRONMENTAL HEAT LOAD TO SUIT" : INPUT 05(8) 

"FAN POWER TO AIRSTREAM" : INPUT 06(13) 

"PUMP POWER HEAT LOAD" : INPUT 06(14) 

"TOTAL PRESSURE - GAS" : INPUT PI (2) 

" C02 PARTIAL PRESSURE INTO SUIT" : IWPUT P3(l) 

"FAN VOLUME FLOW RATE" : INPUT C5 
"PUMP FLOW RATE CH20)": INPUT Wl(l) 

"UA OF LCG (H20)": INPUT U1 

"COOLANT H20 INLET TEMPERATURE TO SUBLIMTOR" : INPUT T4( 


120 PRINT 
130 PRINT 
140 PRINT 
150 PRINT 
160 PRINT 
170 PRINT 
180 PRINT 
190 PRINT 
200 PRINT 
8 ) 

210 PRINT 
220 PRINT 
230 PRINT 
240 PRINT 
250 PRINT 
260 PRINT 
270 PRINT 
280 PRINT 
290 PRINT 
300 PRINT 
310 PRINT 
320 PRINT 
330 PRINT 
340 PRINT 
350 PRINT 
360 PRINT 
P4( 1 ) 

370 PRINT 
380 PRINT 
390 PRINT 
400 PRINT 
410 PRINT 
420 PRINT 


"COOLANT H20 MASS FLOW RATE IN VEHICLE HX": INPUT W2(l) 
"CONDUCTANCE BETWEEN LIOH AND CANISTER" : INPUT Cl 
"CONDUCTANCE BETWEEN GAS AND CANISTER" : INPUT C2 
"CONDUCTANCE BETWEEN CANISTER AND WATER LOOP": INPUT C3 
"CONDUCTANCE BETWEEN LIOH AND GAS STREAM" : INPUT C4 
"MAXIMUM RELATIVE HUMIDITY LEAVING SUIT": INPUT HI 
"GUESS OF FAN INLET TEMPERATURE" : INPUT Tl(13> 

"GUESS OF GAS SUIT INLET TEMPERATURE" : INPUT T1 ( 5) 

"TYPE OF SYSTEM - 1=IVA, 2-EVA" : INPUT K 
"UA OF VEHICLE HX": INPUT U2 

" H20 FLOW TO LIOH CANISTER FOR COOLING" : INPUT Wl(19) 

"POROUS PLATE AREA" : INPUT A1 

"UA OF H20 LOOP IN SUBLIMATOR" : INPUT U3 

"SUBLIMATGR FIN RESISTANCE" : INPUT Rl 

"UA OF GAS LOOP IN SUBLIMATOR" : INPUT U4 

"PRESSURE LOSS BETWEEN SUIT INLET TO FAN INLET" : INPUT 

"HELMET CONDUCTANCE TO INLET GAS STREAM" : INPUT U5 
"INITIAL SKIN TEMP AT REST" : INPUT T3( 1 ) 

"ENVIRONMENTAL HEAT LOAD TO DCM":I!\PUT Q1 (9) 
"ELECTRICAL HEAT LOAD TO DCM" .* INPUT Ql(10) 

"ELECTRICAL HEAT LOAD TO EVCS" : INPUT Q5(2) 
"ENVIRONMENTAL HEAT LOAD TO PLSS" : INPUT Q3(10) 


TABLE II CONT. 


430 PRINT "INC -FACTOR TO ALTER SUBLIMATOR UA (TYPICALLY* 1 ) " .* INPU 
T I 

440 PRINT "DO YOU WISH TO HAVE A DIAGNOSTIC PRINTOUT? (INPUT 0 F 
OR NO OR 1 FOR YES) 11 : INPUT M 
450 SELECT PRINT 215(132) 

460 PRINT "SUITSK2 PROGRAM RUN" : PRINT :PRINT 
470 PRINT SPRINT "PROGRAM INPUT DATA: PRINT 
480 PRINT " QMETT , QENV , OF AN , QPUMP = " ; 0 1 ( 8 ) ! 05 ( 8 > 5 06 < 1 3 > : 06 ( 1 4 > 

490 PRINT " PTOT , PC02, CFM, WH20, UAST=" ; P 1 ( 2) ; P3 ( 1 ) ; C5; W1 (1 ) ; U1 
500 PRINT " TWI , WCPW, CLIOH, CL2=" S T4 (8 ) ; W2 ( 1 > ; Cl S C2 
510 PRINT "CL3,CL4,RH=";C3;C4;Hi 

520 PR I NT " TF I , TSUBO , XKEY , UAX , WLIOH, AP, UAW= " ; T 1 ( 1 3 > *, T 1 ( 5 ) ; K ; U2 ; W 
1 ( 1 9 ) * A 1 " U3 

530 PRINT "RFIN ! ,HAGAS,CDP,UAHEL=":Rl;U4;P4<l);U5 

540 PRINT " TSKINT, GDC MEN , GDC MEL , QEVC S , QPLSEN , i NC * " ,* T3 ( 1 ) ; G 1 ( 9 ) ; Q 

1 ( 10) ; G5(2) ; Q3( 10) ; I 

550 IF M=1 THEN 570 

560 SELECT PRINT 005(64) 

570 02 (8) = 1040* (.0655+. 0252*01(8) MOO) 

580 T3(2)=T3(1>-, 008125*01 (8) 

590 T4(6)~T3(2>"5 

600 V2(8)=. 0001625*01 (S) 

610 V3(8)-1.241*V2(8> :K1=K 
620 R1=R1*I 
630 U4=U4*I 
640 A1=A1*'I 

650 W1 (2) -W1 ( 1 ) : W1 (8) -W1 ( 1 > 

660 T4( 18) =0 
670 Q3 ( 18 ) =0 
680 Q4 ( 8 ) =0 
690 01 (28) =01 (8) 

700 W1 ( 10) =0 
710 Q1 ( 2) -0 
720 05 ( 1 ) =0 
730 K2-K 1 

740 01(2) =Q 1(9) +01(10) 

750 G5 ( 1 ) =Q5 ( 2 ) +G3 ( 10 ) 

760 Q2 ( 7 ) =V3 ( 8 ) *860 

770 02 ( 17)=V3(8>*430 

780 Q3(8) =0 

790 K3=0 

800 S=C 1 +C2+C3 

810 01(7) =Q2( 17 ) +Q2( 8 ) 

820 B=U1/W1(1) 

830 T1 ( 1 ) =T1 (13) +459. 6 
840 T1 (4) =T1 (5) 

850 P2(4)=FN1(T1(4) ) 

860 T 1 ( 4 ) -T 1 ( 4 ) +459 . 6 


8 


TABLE II CONT. 


870 IF P2(4)OPl (2)/2 THEN 890 
880 PS(4)=P1<S)/S 

890 PRINT :PRINT "TA4, PW4=" ; T1 (4) ; P2<4) : PRINT 
900 D1 * ( PI (2) -PE ( 4 ) > *2. 981 /t 1 < 1 ) 

910 VI ( 1 )=D1*C5*60 

920 V3 ( 1 ) =P3( 1 ) / (PI (2) -P2(4) >*44/32*Vl < 1 > 

930 V2( 4 ) =V1 < 1 ) -V3 ( 1 ) 

940 V4 ( 4 ) =C 5*60*P£ ( 4 ) * 1 44/85 . 7/T 1 < 4 ) 

950 VI <4)=V1 (1 ) +V4(4) 

960 FOR N~1 TO SO: PRINT : PR I NT 

970 PRINT "ITERATION CYCLE NUMBER" ; N: PRINT 

980 03 (9) =03(8) 

990 N7-0 

1000 01(7) =02 (17) +02(8) +03(8 ) 

1010 T3 ( 2 ) =T3 ( t ) - . 008 1 25* (01(8) -03 ( 8 ) > 

1020 P2(6)=FN1 <T3(2) ) 

1030 P2(6)=P2(6)*H1 
1040 V2 ( 5 > = V2 ( 4 ) + V2 < 8 ) 

1050 VI ( 5) =V2( 5) +V4<4) +V3( 1 ) 

1060 A2=V4(4)/V2(5) 

1070 R5= ( V4(4 )*85. 7+V3( 1 )*35. 1 +V2< 5 >*48. 3 > /VI < 5 > 

1080 P2<5)=P1(2>*V4(4>/V1<5>*85.7/R5 
1090 PRINT "SUIT INLET CONDITIONS:" 

1100 PRINT ,, PW5,R5,WT5,U45,QSW=";P2(5);R5;V1(5);V4(4);03(8) 

1110 V3(6)=V3(1)+V3(8) 

1120 V2 ( 6 ) - V2 ( 5 ) -V2 ( 8 ) 

1130 IF 01(8) > = 7 49 . 99999 THEN 1680 
1140 IF K3=l THEN 1180 
1150 V4< 6 ) =V4( 4 ) +02 ( 8 ) /1040+Q3<8 ) / 1040 
1160 V 1 < 6 ) = V4 < 6 ) + V2 ( 6 ) + V3 ( 6 ) 

1170 GOTO 1400 
1180 P2( 16) =P2(6) 

1190 FOR Z=1 TO 10 

1200 V4(6) =V1 (6) *P2 ( 16 ) /P 1 ( 2 > *R6/85. 7 
1210 VI ( 16 ) -V4( 6) +V2(6 ) +V3 (6 ) 

1220 IF ABS(V1 ( 16)-V1(6) ) <=V1 <6>*„01 THEN 1240 

1230 GOTO 1250 

1240 IF Z>2 THEN 1310 

1250 VI <6)=V1 ( 16) 

1260 R6=(V4(6)*S5.7+V2(6)*48.3+V3<6)*35. 1)/V1(6) 

1270 NEXT Z 

1280 PRINT "H2Q FLOW SUIT OUT NOT CONVERGING" 

1290 PRINT "WW6,PW6 J WT8=";V4<6>;P2(16);V1(G> 

1300 STOP 

1310 02 ( 16) =1G40*( V4(6) -V4( 4 ) ) 

1320 IF 01 (8) >=749. 99999 THEN 1340 
1330 GOTO 1350 

1340 03 ( 18) =03(8) - <02(16) -02(8) ) ORIGINAL PAGE IS 

OP POOR QUALITY 


TABLE II CONT. 


1350 IF 03(18X0 THEN 1370 
1360 GOTO 1380 
1370 03 ( IS) -0 
1380 Q3(8)=Q2< 16X02(8) 

1330 Q3(9) =03(8) 

1400 03(6 ) =V1 ( 6) *« 22* (T3 (2) -T1 ( 4 ) +459. 6) 

1410 IF 01 (8 >>=749.99999 THEN 1970 
1420 REM LOG HEAT LOAD 

1430 01 (18) ==01 (S> -03(6) +05(8) -03(8) -02(8) 

1440 REM AVERAGE GAS CONSTANT OUT OF SUIT 
1450 R6=(V3(6)*35. 1+V4(6>*85. 7+V2(6)*48. 3) /VI (6) 

146.0 REM SUIT OUTLET PARTIAL PRESSURE 
1470 P3<6)=P1<2)*(V3<6>*35. 1 ) / ( VI <G)*R6> 

1480 P2 ( 1 6 > =P 1 ( 2 ) * ( V4 ( 6 ) *85 » 7 > / < V 1 (6)*R6) 

1490 P4 ( 6 > =P2 ( 6 ) /P2( 16) - 1 
1500 IF P4(6)>. 001 THEN 1520 
1510 GOTO 1530 
1520 IF 03(8X0 THEN 1560 
1530 IF P2 ( 16) >P2 ( 6) THEN 1620 

1540 PRINT 11 WW6, WT6, QSW PW6, PW6M=" ; V4<6) 5 VI (6) 5 03(8) 5 P2( 16) 5 P2( 6 
) 

1550 GOTO 1990 

1560 V4(6)=V4(6)*P2(6)/P2<16> 

1570 03 ( 8 ) = ( V4 ( 6 ) -V4 ( 4 ) ) *1040-02(8) 

1580 01(7) =Q2(8 ) +02 ( 17) +03(8) 

1590 01 ( 18) ==Q1 (8) -03 (6) +05 (8) -03(8) -02(8) 

1600 PRINT ,, WW6,QSW,PP6= ,, ;V4{6);03(S);P4(6) 

1610 GOTO 1160 

1620 V4 ( 6) =V4(6) *P2(6 ) /P2( 16 ) -.001 
1630 03(8) = ( V4(6) ~V4(4) ) *1040-02(8) 

1640 01 (7) =02(8) +02(17) +03(8) 

1650 01 (18) =01 (8) -03(6) +05(8) -03 (8) -02(8) 

1660 PRINT "WW6,QSW= ,, ;V4(6) ”,03(8) 

1670 GOTO 1160 

1680 IF N7> 10 THEN 1970 

1690 V4 ( 6 ) =V4 (4 ) +02 ( 8 ) / 1040+03(8) / 1040 

1700 VI ( 6 ) =V4\( 6 ) +V2 ( 6 ) +V3 ( 6 ) 

1710 03(6) =V1 (6)*. 22* ( T3( 2) -TK 4 ) +459.6 ) 

1720 REM DETERMINE QLCGMET FROM GRAPHICAL DATA 

1730 IF 01 (8X=1 100 THEN 1750 

1740 GOTO 1760 

1750 01 (17) =.65*01 (8) +47. 5 

1760 IF 01 (8) >1100 THEN 1780 

1770 GOTO 1790 

1780 01 (17) =.8438*01 (8) -165.7 

1790 REM DETERMINE TOTAL SWEATING RATE 

1800 03(8) =01 (8) -01 (17) -03(6) -02(8) 

1810 IF 03(8X0 THEN 1830 


ORIGINAL PAGE IS 
OF poor quality; 


10 


TABLE II CONT. 


1880 GOTO 1340 

1830 Ql(17) =01(17) +03(8 ) 

1840 IF 03(8X0 THEN I860 
1850 GOTO 1870 
1860 03(8) =0 

1870 IF ABS(Q3(8) -03(9) )<=. 1 THEN 1910 
1880 03 (9) =03(8) 

1890 GOTO 1680 

1900 REM DETERMINE ACTIVE SWEATING AND INSENSIBLE SWEATING 
1910 R6=(V3(6)*35. 1+V4(6>*85. 7+V2(6)*48. 3) /VI (6) 

1980 P2( 1S)=P1 (2)*V4(6)*85. 7/VI <6)/R6 
1930 IF P2( 16)<=P2(6) THEN 1970 
1940 N7=N7+1 

1950 PRINT "QSW,QSWA,PW6,PW6M,N7= ,, ;G3(8>;Q3(18);P2(16);P2(6);N7 

1960 GOTO 1180 

1970 01(18) =01(17) +QS(8) 

1980 P3(6)=P1(2) *V3 ( 6 ) #35 . 1/V1(6)/R6 

1990 T1 ( 6 ) =T3 ( 2 ) +459. 6 

2000 P2( 15) =P1 (2)*V4(6)*85. 7/ VI (6)/R6 

2010 PRINT " WW6, WT6, QSW, PW6N, PW6M= " ; V4 ( 6 ) ; VI ( 6 ) ; 03 ( 8 ) ; P2 ( 1 5 ) ; P2 ( 
6) 

2020 REM LOOP ON LI OH OUTLET TEMPERATURE 

2030 IF N=1 THEN 2050 

2040 GOTO 2060 

2050 T4(2) =T1 (6) -469. 6 

2060 T3( 3) =T4(2) +459. 6 

2070 V3 ( 16 ) =V1 (6) 22 

2080 Y( 1 , 1 XC1+C4 

2090 Y(2, 1 ) =C4 

2100 Y(3, 1)=C1 

2110 Y(4, 1 )=0 

2120 Y( 1 , 2) =-C4 

2130 Y (2, 2) =-C4~V3 ( 16) 

2140 Y( 3, 2 ) =0 
2150 Y(4, 2) =V3( 16) 

2160 Y( 1 , 3) =-Cl 

2170 Y(2, 3) =0 

2180 Y( 3, 3 ) =-S 

2190 Y(4, 3) =C2 

2200 Y ( 1 , 4 ) =0 

2210 Y ( 2 , 4 ) =0 

2220 Y( 3, 4 ) =02 

2230 Y(4, 4 ) =-V3( 16 ) -C2 

2240 Z(l, 1)=Q2(7) 

2250 Z(2, 1)=-V3(16)*T1(6) 

2260 Z ( 3, 1 ) =-C3*T3 ( 3 ) 

2270 Z (4, 1 ) =0 
2280 MAT I=INV(Y),M 
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£290 IF M~Q THEN 2310 
2300 GOTO 2330 

2310 PRINT "LIOH TEMPS NOT CONVERGING" 

2320 STOP 
2330 MAT X=I*Z 
2340 T5( 1 ) =X ( 1, 1 ) 

2350 T5(2) =X (2, 1 > 

2360 T6(1)*X (3,1) 

2370 T1 (7 ) =X (4, 1 ) 

2380 PRINT "LIOH TEMPS CONVERGED" 

2390 PRINT "QLS,WT6,TA7,TA6,TWP0=";Q£(7);V1(6);T1(7);T1(6);T4(2) 

£400 PRINT m TLI,TL0,TC,TA7=";T5<1);T5(2)*,TG(1>;T1(7) 

2410 REM FAN INLET CONDITIONS 
2420 VI ( 7 ) =V1 (6) -V3(8) +02 ( 17 ) / 1060 
2430 V3 ( 7 ) =V3 < 6 ) - V3 ( 8 ) 

2440 IF V3(7) <0 THEN 2460 
2450 GOTO 2470 
2460 V3(7) =0 

2470 V4 (7) =V4< 6) +02 ( 17) / 1060 
2480 V2(7 ) =V2(6 ) 

2490 R7= ( V2 ( 7) *48. 3+V3 (7 ) *35. 1+V4(7>*85. 7) /VI (7) 

2500 Pl<l)=Pl<2>-P4(i>*<Vl<7m.30)/Pl<£>*3.8/560*Tl<7> 

2510 D7=P1 ( 1 )*144/R7/T1 (7) 

2520 P3(7)=V3(7)/V1 (7)*35. 1/R7*P1 (2) 

2530 P2 ( 7 ) =V4 (7 ) /VI ( 7 ) *85 u 7/R7*P 1 < 2) 

2540 REM CALCULATED FAN FLOW RATE 
2550 V1(17)=D7*C5*60 

2560 PRINT 41 WC7 , WW7 , W07 , R7, P IN= " ; V3 ( 7 ) ; V4 ( 7 ) 5 V2 ( 7 ) ; R7; P 1 ( 1 ) 

2570 PRINT "RH07,PC7,PW7,WT7,WT7C=";D7;P3<7>;P2(7);V1(7);V1(17) 

2580 REM FAN OUTL.ET TEMP 

2590 T1(2)=T1(7)+QG(13>/(V1 (17)*. 22) 

2600 REM LIOH SENS HEAT LOAD TO GAS 
2610 Q3(19>=V1<6>*.22*(T1(7>~T1<6) ) 

2620 REM TOTAL GAS HEAT LOAD SENS 
2630 G3( 3 ) -Q3< 19 > +Q6 (13) +Q3(6) 

2640 IF Q3 ( 3) <=0 THEN 2660 

2650 GOTO 2670 

2660 Q3( 3) =03(19) +06(13) 

2670 01 (7 ) "02(8 ) +Q3 ( 8 ) +02 ( 17 ) 

2680 01 (3) =03(3) +Q1 (7) 

2690 REM WATER LOOP SUIT INLET TEMP 
2700 T4 ( 16 > =T3 ( 2 ) -01 ( 18 ) /W1 ( 1 ) / < 1 -EXP ( -B ) ) 

2710 PRINT "TWSI,QLCG,PIN=";T4(16):Q1(18);P1(1) 

2720 REM SUIT OUTLET TEMP-H20 
2730 T4(6) =T4< 16) +01 ( 18) /W1 ( 1 ) 

2740 T4(2) =T4(6) +(06(14) +05 ( 1 ) +02(7) -03 ( 19) ) /W1 ( 1 ) 

2750 T4( 19) =T4(2)+ (02(7) -03(19) )/Wl(19> 

2760 REM TOTAL H20 LOOP HEAT LOAD 
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2770 Q1(23)=Q1( L8) +Q6< 14) +Q2( 7 > -Q3< 19) +05 ( 1 ) +01 (2) 

2780 REM SUBLIMATOR WATER INLET TEMP 
2790 T4( 13>=T4( 16) +01 <23>/Wl ( 1 ) 

2800 PR I NT " TWS I , TWSO , TWSB I , TWPO , QWT= " ; T4 ( 16 ) T4 C 6 ) T4 ( 1 3 ) *, T4 < 2 ) 
;G1<23) 

2810 01(13) =01 ( 23 ) +01 ( 3 ) 

2820 IF K2=2 THEN 4120 
2830 REM TOTAL HEAT LOAD 
2840 01 (13) =01 (23) +01 (3) 

2850 REM USING WATER LOOP COOLING 
2860 W1(S)=W1(1> 

2870 T2=0 

2880 FOR N3=l TO 10 

2890 C=U2/W1 (8»( 1-W1 <8)/W2( 1 ) ) 

2900 El=< l-EXP(-C) )/(l-Wl(8)/W2(l)*EXP(-C> ) 

2910 T4( 18)=T4( 13>-E1*(T4< 13)-T4(8) ) 

2920 G1(53)=W1(8)*(T4(13)-T4(18) ) 

2930 PRINT "QT , QTC , TWXO, EFX , WHX = " ; 01 ( 13) 01 ( 53) ; T4( 18) ; El *, W1 ( 8) 
2940 REM IS THE 0 IN TOLERANCE? 

2930 IF 01 (13) >01 (53) THEN 3730 
2960 IF 01 (13) >01 (53) THEN 3080 
2970 IF N3=l THEN 2990 
2980 GOTO 3000 
2990 01 (52) -01 (53) 

3000 IF ( ABS(Q1 (53) -01(13) )- .01*01 (13) > <=0 THEN 3080 
3010 W1 (8)«W1(8)*Q1(13)/Qi(53)*.75+.25*W1 (8) 

3020 NEXT N3 

3030 REM LOOP NOT CONOIRGED 
3040 PRINT "HX LOOP NOT CONVG" 

3050 PRINT i, QT,QTC = ,, ;Q1(13);Q1(53) 

3060 STOP 

3070 REM CAN HX DO JOB? 

3080 IF W1 (8)>W1 ( 1 ) THEN 3100 
3090 GOTO 3130 

3100 PRINT " HX NEEDS TOO MUCH FLOW" 

3110 PRINT ,1 WW,WHX = M ;W1(1);W1<8) 

3120 STOP 

3130 REM SUBLIMATOR PERFORMANCE- GAS TEMPS 

3140 REM MASS FLOW RATIO 

3150 W4=V1(17)*G1 <3)/Q3(3)/Wl <8)*.22 

3160 T2= (T4( 13) +T4( 18) ) /2 

3170 IF T2>S5 THEN 3220 

3180 IF T2<>55 THEN 3200 

3190 GOTO 3210 

3200 IF T2>47. 5 THEN 3270 

3210 IF T2<=47.5 THEN 3320 

3220 U3= (2. 5785*W 1 ( 8 ) + 1 24 . 01 ) * I 

3230 IF W1 (8X144.63 THEN 3250 
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3240 GOTO 3260 

3250 □3=419. 52# I 

3260 GOTO 3370 

3270 U3=(2. 5*W1(8)+98)*I 

3280 IF WK8K124.57 THEN 3300 

3290 GOTO 3310 

3300 U3=41 0.69*1 

3310 GOTO 3370 

3320 U3= (2. 212*W1 (8) +121 . 08 )*I 
3330 IF W1 (8X128 THEN 3350 
3340 GOTO 3360 
3350 U3=406*I 
3360 GOTO 3370 

3370 U9=U4*Q1 ( 3 ) /Q3 ( 3 ) *U3/ ( U4*Q1 ( 3 ) /Q3 ( 3 > +U3 ) 

3380 D=U9/(W4*W1(8) )*<1~W4) 

3390 E2= ( 1 -EXP ( -D ) > / < 1 ~W4*EXP ( -D ) ) 

3400 IF W40.125 THEN 3420 
3410 GOTO 3430 
3420 E2= « 98*E2 

3430 T1 < 22) =T4( 18 ) +Q3( 3) / V 1(17)/. 22/E2+459. €■ 

3440 T1 (24) =T1 (22) -Q3( 3) /VI ( 17 > /. 22 
3450 PRINT "TAVG,UAW,WHX=";T2;U3;W1(8) 

3460 PRINT 11 TWXO, TA2C , TA4C , EFSX , TA2= " ; T4 ( 18 ) ; T1 ( 22 ) ; T 1 ( 24 ) 5 E2; T1 
(2) 

3470 PRINT ,, WR,QSSUB,D,C = ";W4;Q3(3);D;C 
3480 PRINT 

3490 T4 ( 14) =T4( 18) +01 ( 3) /W1 (8) 

3500 W1 < 12 ) =W1 ( 1 ) -W1 (8 ) 

3510 IF W1 ( 12) >0 THEN 3530 

3520 GOTO 3540 

3530 IF K3= 1 THEN 3550 

3540 GOTO 3680 

3550 Q3( 18) =0 

3560 W1(S)=W1(1) 

3570 W1 ( 12 ) =0 
3580 Q1(6)=G1(8) 

3590 01(8) =Q1 (6 ) *Q1 ( 52) /Q1 ( 1 3) 

3600 IF ABS(Q.l (6X01(8) )> = 5 THEN 3620 

3610 GOTO 3630 

3620 01 (8) =01 (6) +5 

3630 IF 01 (8) >01 (28) THEN 3650 

3640 GOTO 3660 

3650 01 (8) =01 (28) 

3660 K3=0 
3670 N=1 
3680 01(33) =0 
3690 01 (23) =01 (3) 

3700 PRINT ,, QMTI,QMT=“;Q1(6);Q1(8) 
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3710 REM GQ TO CONVERGENCE CHECK 

37S0 GOTO 5340 

3730 IF K3>0 THEN 4010 

3740 REM FIRST ITERATION SET KEY 

3750 IF N>3 THEN 3770 

3760 GOTO 3790 

3770 K3=l 

3780 REM SUIT OUTLET H20 PRESS 
3790 P2( 16) =P2(6)*H1 
3800 REM LOOP TO FIND H20 FLOW RATE 
3810 FOR N4=l TO 10 

3820 V4 < 6 ) =P2 ( 1 6 > /P 1 < 2 ) *V1 ( 6 > *R6/85 . 7 
3830 V 1 < 16 ) = V4 £ 6 > + V2 ( 6 ) +V3 ( 6 ) 

3840 IF ( ABS( VI < 16) -VI (6 ) ) -• Q1*V1 (6) ) <=0 THEN 3940 
3850 REM NOT CONVG NEXT ITER. 

3860 VI (6) =V1 (16) 

3870 R6=(V4(6>*85.7+V2<6)*48. 3+V3(6)*35. 1 )/Vl (6) 

3880 NEXT N4 

3890 REM LOOP NOT CONVG 

3900 PRINT "H20 SUIT OUTLET NOT CONVG" 

3910 PRINT "WW6,PW6,WT6=" ; V4(6> ;P2< 16) , VI (6) 

3920 STOP 

3930 REM SWEAT RATE 

3940 02(16) *1040* ( V4 C 6 ) -V4 ( 4 ) ) 

3950 03(8) “02 (16) -02 ( 8 ) 

3960 IF 03(8X0 THEN 3980 
3970 GOTO 3990 
3980 03(8) =0 
3990 GOTO 5420 

4000 REM REDUCE METABOLIC HEAT LOAD 

4010 Q1(8)=Q1<8)-S 

4020 V2(8) -.0001625*01(8) 

4030 02 (8 >*1040* (.0655+. 0252*01 (8) MOO) 

4040 T3 < 2 > =T3 ( 1 ) - . 008125*01 ( 8 ) 

4050 V3 ( 8 ) = 1 . 24 1 *V2 ( 8 ) 

4060 02(7) * V3 ( 8 ) *860 
4070 02 ( 17 ) =V3 (8 ) *430 
4080 03 ( 18) =0 

4090 PRINT " QMT REDUCED: GMT, QML= " ; 01 (8 ) ; 02 ( 8 ) '■ 
4100 GOTO 5420 

4110 REM LOOP ON SUBLIMATOR PERFORMANCE 
4120 W1 (8) "W1 ( 1 ) 

4130 REM LOOP ON FLOW RATE 
4140 FOR L2=l TO 20 
4150 T3=0:N9=0:T2=0 
4160 T4(3)=Q 1(23) /W 1(8) 

4170 REM MASS FLOW 

4180 W1 (3>=(Q1 (23) +01 (3) )/T4(3) 
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4190 REM Q/A PLATE 

4500 Q1 (43) - (Q1 (23)+Ql (3) )/Al 

4510 G 1 = 3. 818E --07*01 (43) +.001585 

4550 U6 = 3 « 96* A 1 / G 1 

4530 REM SUBLIMATOR WATER UA 

4540 07=03*06/ ( U3+U6 ) 

4550 N9-N9+ 1 : IF N9-10 THEN 4600 

4560 PRINT " TOLD, TAVG , NNN= " ; T3; T5? N9 

4570 PRINT " UASW , UAW , UAG , WMCP = " ; 07 ; U3 ; U6 ; W 1 ( 3 ) 

4580 PRINT "DTW, GAP, GAP= " ; T4 ( 3) ; 01 (43) ; G1 
4590 REM OUTLET TEMP 

4300 T4 ( 14 ) =35+ ( T4 ( 13 ) -35) /EXP ( U7/W1 ( 3 ) ) 

4310 IF T4( 14) <33 THEN 4330 
4350 GOTO 4340 
4330 T4< 14) -33 

4340 IF ABS ( T3-T4 ( 1 4 > ) < = . 1 THEN 4600 

4350 T3=T4( .14) : T5-- (T4( 13) +T4( 14) ) /5 

4360 IF T5<=75 THEN 4370: GOTO 4380 

4370 IF T5>65 THEN 4430 

4380 IF T2<=65 THEN 4390: GOTO 4400 

4390 IF T5>55 THEN 4470 

4400 IF T2<=55 THEN 4410: GOTO 4450 

4410 IF TS>47. 5 THEN 4510 

4450 IF T2<=47.5 THEN 4550 

4430 03= (5. 75*W1 ( 8 ) +138 ) *1 

4440 IF ( U 3-468* I ) <0 THEN 4450: GOTO 4460 

4450 U3=468*I 

4460 GOTO 4540 

4470 U3=<2.6*Wi(8)+118)*I 

4480 IF (03-450*1 )<0 THEN 4490: GOTO 4500 

4490 U3=420*I 

4500 GOTO 4240 

4510 U3=Ce.425*Wl (8>+110>*I 

4520 IF (U3-410*IX0 THEN 4530: GOTO 4540 

4530 U3=410*I 

4540 GOTO 4240 

4550 03= (2. 121*W1(8)+121„08)*I 

4560 IF (03-406*1X0 THEN 4570: GOTO 4580 

4570 03=406*1 

4580 GOTO 4240 

4590 REM CALCULATED HEAT LOAD 
4600 Q1 (53XW1 (3)*(T4( 13)-T4( 14) ) 

4610 IF L2=i THEN 4620: GOTO 4630 
4620 01(52) =01 (53) 

4630 REM IS 0 WITHIN TOLERANCE ? 

4640 PRINT H TWSIO,TWSBI,QTC,WHX=" ;T4( 14);T4(13) ;Q1(53);W1 (8) 
4650 IF 01 (13) >01 (53) THEN 4770 
4660 IF 01 (13) >01 (53) THEN 4680 
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4670 GOTO 4690 
4680 09= (OK 13) -Q1 (53) )/4 
4690 03 ( 8 ) =03 < 8 ) +09 
4700 09=0 

4710 IF 01 (13) >01(53) THEN 4750: GOTO 4740 
4730 IF N=3 THEN 4730: GOTO 4740 
4730 K3=l 

4740 IF 01 ( 13) >01 ( 53) THEN 4910 

4750 IF ( ABS ( 01 ( 53 ) -01 ( 13) ) - . 01*01 (13) ) <=0 THEN 4910 

4760 GOTO 4S30 

4770 IF L3>1 THEN 4910 

4780 GOTO 4660 

4790 IF N=1 THEN 4810 

4800 GOTO 4740 

4810 K3= 1 

4830 GOTO 4740 

4830 REM NOT IN TOL- NEW VALUE OF Wl<8) 

4840 W1 (8)=W1 (8)*<G1 ( 13)/Q1(53) ) 

4850 NEXT L3 

4860 REM LOOP NOT CONVG- PRINT 
4870 PRINT "HX LOOP NOT CONVG " 

4880 PRINT "QT„QTC = " ’,01(13) ;Q1 (53) 

4890 STOP 

4900 REM IF FLOW TOO HIGH 
4910 IF W1 (8) >W1 ( 1 ) THEN 4930 
4950 GOTO 4960 

4930 PRINT 11 HX NEEDS TOO MUCH FLOW" 

4940 PRINT "WW„WHX = "’,W1C1);W1(8) 

4950 STOP 

4960 W1 ( 15) =W1 ( 1 ) -W1 (8 ) 

4970 IF W1 ( 13) >0 THEN 4980: GOTO 4990 
4980 IF K3=l THEN 5000 
4990 GOTO 5050 

5000 09=0: Q3< 18)=0:W1 (8)=W1 (1 ) :W1 (13) =0:01 (6)=Q1 (8) 

50 10 01(8) =01 ( 6 ) *01 ( 53 ) /G1 ( 1 3 ) 

5030 IF ABS (01 ( 6 ) -01 ( 6 ) *0 1 ( 53 ) /Q 1 ( 13) ) > = 5 THEN 5030:G0T0 5040 
5030 01 (8) =01 (6) +5 
5040 N-l 

5050 W1 ( 10) =01 f 13)/ 1037:01 ( 33 ) =Q1 ( 13) 

5060 PRINT " WB VPS , WFEED , QTSBL , GMT , QMT I = " ; W 1 ( 1 3 ) ; W 1 ( 1 0 ) ; 0 1 ( 33 ) Q 1 
(8) : Q) (6) 

5070 REM A IRS IDE OF SUBLIMATOR - WALL TEMP - PRIMARY 
5080 T4 ( 39 ) =33+ CT4 ( 1 3 ) -T4 ( 14 ) ) /LOG ( ( T4 ( 1 3 ) -33 ) / ( T4 ( 14 ) -33 ) ) 

5090 REM METAL TEMP 

5100 T6 (3) =33+01 ( 13)/U6+Q1 ( 13) /Rl+3. 63- 1 . 587E-03*Q1 (13) 

5110 T 1 < 1 4 ) =T 1(4) 

5130 R=LOG( 10*P3(7) ) : T3(7) =FN3(R ) +459. 6 
5130 U 1 ( 4 ) =U4*Q 1 < 3 ) /Q3 ( 3 ) 
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5140 FOR N5=l TO 10 

5150 W3(1)=Q1(7)/(T2(7)-T1(14) > 

5160 T1 ( 24 > =T6( 2) +01 ( 7 ) /W3 ( 1 ) / (EXP ( U1 (4 ) /W3C 1 ) ) -1 ) +459. 6 
5170 PRINT " TA4A,TA4C=" ;T1 < 14) ;T1 (24) 

5180 IF ABS(T 1 ( 14)-T1 (24) ><=, 1 THEN 5260 
5190 REM GAS OUTLET TEMP 
5200 T1 ( 14)=T1 (24) 

5210 NEXT N5 

5220 PRINT "SUBL GAS TEMP NOT CQNVG" 

5230 PRINT " TD7 , WM , TA4C , TM , TWL.= " ; T2(7) J W3( 1 ) ; T1 (24) 5T6(2) ; T4(29) 
5240 STOP 

5250 REM GAS INLET TEMP 

5260 T1 (22) =T1 (24) +Q3( 3) /(VI ( 17)*. 22) 

5270 PRINT " TA2C, TA2=" 5 T1 ( 22) ; T1 (2) 

5280 PRINT" TD7 , WM TA4C , TM , TWL = " ; T2 ( 7 ) ; W3 ( 1 ) ; T 1 ( 24 ) ; T6 ( 2 ) ; T4 ( 29 ) 
5290 T4( 18) =T4( 13) 

5300 IF 01(13) >Q1 ( 53 ) THEN 5320 

5310 GOTO 5330 

5320 IF K3=l THEN 4010 

5330 REM IS THIS THE FIRST ITERATION ? 

5340 IF N<=4 THEN 5420 

5350 02 ( 3 ) = ( V4 ( 7 ) -V4 ( 4 ) )*1040 

5360 PRINT " QLC , QLT- " ; 02 ( 3 ) I 0 1 ( 7 ) 

5370 REM NO- IS SYSTEM CONVERGED ? 

5380 IF ( ABS ( T 1 ( 24 ) -T 1 ( 4 ) ) - . 2 ) >0 THEN 5420 
5390 T3( 12)=T3( 1 >-» 008125* (01 (8)-Q3<8) ) 

5400 IF (ABS ( T3 ( 12) ~T3 (2 ) ) -1 ) <=0 THEN 5640 
5410 REM NOT CONVERGED, NEW ITERATION SET UP 
5420 Vl(.l) =V1 ( 17 ) -V4 ( 7 ) f-V4 ( 4 ) 

5430 T1(4)=T1(24)*.5+.5*T1<4> 

5440 IF T1 (4X491.6 THEN 5450: GOTO 5460 
5450 Tl(4)=491.6 
5460 T1 ( 1 )=T1 (7) 

5470 REM DEW POINT LEAVING HX 
5480 P2(4) =FN1 (T1 (4) -459. 6) 

5490 IF P2(4)>Pl(2)/2 THEN 5510 

5500 GOTO 5520 

5510 P2 (4) =P1 ( 2 ) /2 

5520 REM H20 VAPOR LEAVING HX 

5530 R4=(V2(7)*48.3+V3(7)*35. 1+V4<4)*85. 7) /VI ( 1 ) 

5540 V4(4)=V1 (1 )*P2(4)/P1 (2)*R4/85. 7 
5550 REM C02 LEAVING HX 
5560 V3( 1 ) -P3( 1 ) /PI (2)*V1 ( 1 )*R4/35. 1 
5570 V2(4)=V1 ( 1 )-V4(4) -V3( 1 ) 

5580 PRINT ,! WG1,TA4,TA1,PPW4,R4-";V1(1);T1(4);T1(1);P2(4);R4 
5590 PRINT ,, WW4,WCi,W01--";V4(4);V3(l);V2C4) 

5600 NEXT N 
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Si 

! 5610 REM TEMP LOOP NOT CONVERGED - PRINT 

j 5620 PRINT "LOOP MOT CONVERGED" 

j* 5630 STOP 

j 5640 PRINT "LOOPS CONVERGED. LOAD SECOND PROGRAM, SUITOUT, BY PU 

1 TTING THE APPROPRIATE TAPE CARTRIDGE INTO THE TAPE DRIVE AND P 

RESTING THE CONTINUE/RETURN KEYS." 
i 5650 PRINT 

] 5660 REM GO TO "SUITOUT" FOR DATA OUTPUT 

5670 STOP 

5680 LOAD "SUITOUT" 



TABLE III 


SUITSK2 AND SUITOUT BASIC TO FORTRAN 
VARIABLE CROSS REFERENCE 


BASIC 

FORTRAN 

BASIC 

FORTRAN 

A1 

AP 

P2 (7) 

PW7 

A2 

AHS1 

P2 (15) 

PW6N 

B 

B 

P2 (16) 

PW6 

C 

C 

P3(l) 

PCI 

Cl 

CL1 

P3(2) 

PC2 

C2 

CL2 

P3(6) 

PC6 

C3 

CL 3 

P3(7) 

PC7 

CA 

CLA 

P4(l) 

CDP 

C5 

CFM 

P4(6) 

PP6 

C6 

CP 6 

Q6 

QMETTW 

D 

D 

Q8 

QLATST 

Dl 

RHOl 

Q9 

DQSW 

D7 

RH07 

Ql(2) 

QDCM 

El 

EFX 

Ql(3) 

QSTOT 

E2 

EFSX 

Ql(6) 

QMTI 

G1 

GAP 

Ql(7) 

QLT 

HI 

RH 

QK8) 

QMT 

I 

INC 

Ql(9) 

QDCMEN 

K 

RKEY 

Ql(10) 

QDCMEL 

Kl 

XKEY 

Ql(13) 

QT 

K2 

IKEY 

Ql(17) 

QLCGM 

K3 

KEY 

Ql(18) 

QLCG 

N9 

NNN 

Q1 (23) 

QWT 

Pl(l) 

PIN 

Ql(28) 

QMTO 

Pl(2) 

PT 

Ql(33) 

QTSBL 

P2 (4) 

PWA 

Q1(A3) 

QAP 

P2(5) 

PW5 

Q1 (52) 

QTCI 

P2(6) 

PW6M 

Ql(53) 

QTC 


TABLE III CONT. 


BASIC 

FORTRAN 

BASIC 

FORTRAN 

Q2 (3) 

QLC 

Tl(6) 

TA6 

Q2 (7) 

QLS 

Tl(7) 

TA7 

Q2 (8) 

QML 

T1 (13) 

TF1 

Q2 (16) 

QLMAX 

Tl(14) 

TA4A 

Q2 (17) 

QLL 

T1 (22) 

TA2C 

Q2 (18) 

QLSTW 

T1 (24) 

TA4C 

Q3 (3) 

QSSUB 

T2 (5) 

TD5 

Q3(6) 

QS6 

T2 (6) 

TD6 

Q3(8) 

QSW 

T2(7) 

TD7 

Q3 (9) 

QSWO 

T3 (1) 

T SKINT 

Q3(10) 

QPLSEN 

T3 (2) 

TSK 

Q3 (18) 

QSWA 

T3 (3) 

TSO 

Q3 (19) 

QLSG 

T3 (12) 

TSKC 

Q4(8) 

QSTR 

T4(2> 

TWPO 

Q5 (1) 

QEV2 

T4(3) 

DTW 

05(2) 

QEVCS 

T4(6) 

TWSO 

Q5(8) 

QENV 

T4(8) 

TWI 

Q6(13) 

QFAN 

T4 (13) 

TWSBI 

Q6(14) 

QP 

T4 (14) 

TWSBO 

R1 

RFIN 

T4 (16) 

TWSI 

R4 

R4 

T4(18) 

TWXO 

R5 

R5 

T4 (19) 

TWLO 

R6 

R6 

T4(29) 

TWL 

R7 

R7 

T5(l) 

TLI 

S 

SU 

T5 (2) 

TLO 

T1 

THEL 

T6(l) 

TC 

T2 

TAVG 

T6 (2) 

TM 

T3 

TOLD 

U1 

UAST 

Tl(l) 

TAI 

U2 

UAX 

T1 (2) 

TA2 

U3 

UAW 

Tl(4) 

TA4 

U4 

UAGAS 

Tl(5) 

TSTI 

U5 

UAHEL 
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TABLE III CONCLUDED 


BASIC 

FORTRAN 

BASIC 

U6 

UAG 

Wl (2) 

U7 

UASW 

Wl(3) 

U8 

UASC 

Wl (8) 

U9 

UASX 

Wl(10) 

Ul(4) 

UAGS 

Wl(12) 

VI 

V02A 

Wl(19) 

V2 

V SLURP 

W2 (1) 

V4 

V4 

W3(l) 

V5 

V5 


Vl(l) 

WG1 


VI (4) 

WT4 


VI (5) 

WT5 


Vl (6) 

WT6 


VI (7) 

WT7 


VI (16) 

WT6C 


Vl(17) 

WT7C 


V2(4) 

WOl 


V2(5) 

W05 


V2(6) 

W06 


V2(7) 

WO 7 


V2 (8) 

WO 2 A 


V3(l) 

WC1 


V3(6) 

WC6 


V3(7) 

WC7 


V3 (8) 

WC02 


V3 (16) 

WCP 


V4 (4) 

WW4 


V4 (6) 

WW6 


V4(7) 

WW7 


W4 

WR 


W5 

WSEP 


Wl(l) 

WW 



FORTRAN 

WS 

WMCP 

WHX 

WFEED 

WBYPS 

WLIOH 

WCPW 

WM 


SUITOUT COMPUTER PROGRAM 


Pile Name "SUITOUT" 

Abstract "SUITOUT" provide.s the printout for the data calculated 
in the "SUITSK2" program. Figures 2 and 3 show tabu- 
larized data output and the flowchart data format pro- 
vided by the program. This program is the second half 
of the Hamilton Standard "SUITSK2" program which has 
been converted from Fortran to Basic for use on the 
Wang 2200/Flatbed plotter system. 

Program Description 

The program is automatically loaded by the "SUITSK2" program. The 
user is required only to load the program tape cartridge into the tape 
drive unit as instructed by the "SUITSK2" program. The program pro- 
vides data output for the variables calculated in the "SUITSK2" program 
by printing the tabularized output as exemplified in Figure 2. Table 
IV identifies the output data definitions and units. In addition, an 
optional flowchart is provided as shown in Figure 3. If the user 
desires the flowchart, the CRT instructs him to load a blank flowchart 
(Figure 4) onto the flatbed plotter. The plotter then plots the data 
onto the blank flowchart. Both of the examples of output data shown 
in Figures 2 and 3 are the result of the input data shown in Figure 1. 

Because the Wang 2200 system does not automatically round off numbers, 
a routine is included in the program which performs that function. 
Therefore, the data provided in both the tabularized format and the 
flowchart are rounded off values. Also, a special function, DEFFN2, 
is included in the program to calculate a saturation temperature in 
degrees F based upon the log of 10 times the saturation pressure in 
psia. The function is based upon a subroutine "KANDK" used by the 
original Fortran version of the program. A listing of program "SUITOUT" 
is provided in Table V with a cross reference between the original 
Fortran variable names and the Basic variable names provided in Table 
III. 
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FIGURE 2 SUITOUT PROGRAM EXAMPLE 
TABULARIZED OUTPUT 


PROGRAM OUTPUT DATA” 

suit conditions; 

WTQTI j, WH501 , WCG2I , PH20I ? PC02I , TGASI 3 TDPI^ 31.67 -54 

3.00000000E-02 -1995 .5005 53.5 53.1 

WTOTO, WH200, WCG2Q, PH500, PC050, TGASO, TDPO= 31.83 .36 -53 

.3063 4.0894 87.1 65.1 

THEL. , TSK IN , TLCG I , TL_CGO= 55.1 87.155 74.85 78.55 

WSUITj W05, VQ2, WC05, VSI = 185 .163 3. 06000000E-02 .5016655 

6-01 

QMANS, QMANRL, QSW, QMANTL, QLCG , QLCGT= 537.4 133.64 0 133.6 

658.9 628.9 

LI OH CONDITIONS: 

WTOTO, WH2D0, WC050, PH20Q, PC020, TGASO, TDPO= 31.71 .45 

3. OOOOOOOOE-02 .3752 .4944 95.9 71 

TLIOH, TGO, TCASEn TH20I , TH200, PTOTO= 101.5 99.8 87.4 79.9 

85.5 15.0798 

QL.S , QLL , QLSG , QLH20= 173.4 86.7 61.2 115.2 


8UBLIMAT0R CONDITIONS; 


WSUBL, WBYPS, WFEED, WSEP , VSLURP, VSUBQ= 42. 54 142.46 

7 . 00000000E-02 5 . 97 

TG AS I , TH20 1 , TH200 TVHX I , TVHXO= 97 .4 80.01 57 „ 46 

QSIJBL, QSUBS, QSUBVT, QWTOT , QTSUB , QVHXT- 520.4 309. 7 

1484'. 1 

REVISED PROGRAM INPUT DATA: 

GsMETT, QEN V OF AN ,, QPUMP= 1000 0 11 145 

PTOT, PC02, CFM ? WH20, UAST= 15.2 9. 67000000E-03 6.6 

TWI,WCPW,-CLI0H,CL2= 45 600 6 3.2 

CL3i CL4, RH- 15 50 .9 

TFI,TSUBO,XKEY,UAX,WLIOH,AP !( UAW= 115 SO 1 864.6 

419. 52 


0 .2103 

45 45 

530 530 


185 -60 
21 .2361 


RFIN, HAG AS, CDP, UAHEL- 471.6 16.28 5. 45000000E-03 .4 

TSK I NT, QDCMEN, QDCMEL,, QEVCS,, QPLSEN, INC- 95. 25 0 27 41 0 1 


0 




OR!GIN AL PAGE IS 


SHUTTLE EMU FLOW CHART 

CONDITIONS i MODE- IV 

METABOLIC RATE- 1000 BTU/HR 
ENVIRONMENT- 0000 BTU/HR 


T= 53.3 
DF= 53.3 
V- 5.37 
M= 31.50 

a 



FIGURE 3 SUITOUT PROGRAM EXAMPLE FLOWCHART OUTPUT 
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NAME 

WTOTI 

WH20I 

WC02I 

PH20I 

PC02I 

TGASI 

TDPI 

WTOTO 

WH200 

WC020 

PH200 

PC020 

TGASO 

TDPO 

THEL 

TSKIN 

TLCGI 

TLCGO 

WSUIT 

W02 

V02 

WC02 

VSI 

QMANS 

QMANLR 

QSW 

QMANTL 

QLCG 

QLCGT 

TLIOH 

TGO 

TCASE 

TH20I 

TH200 

PTOTO 

QLS 

QLL 

QLSG 

QLH20 LI OH 

WSUBL 

WBYPS 

WFEED 

WSEP 

V SLURP 

VSUBO 

TVHXI 

TVHXO 

QSUBL 

QSUBS 

QSUBT 

QWTOT 

QTSUB 

QUHXT 


TABLE IV SUITOUT OUTPUT DATA DESCRIPTION 


DESCRIPTION 


TOTAL GAS FLOW RATE AT INLET 
WATER VAPOR FLOW RATE AT INLET 
C02 FLOW RATE AT INLET 
PARTIAL PRESSURE OF H20 AT INLET 
PARTIAL PRESSURE OF C02 AT INLET 
GAS TEMPERATURE AT INLET 
DEW POINT TEMPERATURE AT INLET 
TOTAL GAS FLOW RATE AT OUTLET 
WATER VAPOR FLOW RATE AT OUTLET 
C02 FLOW RATE AT OUTLET 
PARTIAL PRESSURE OF H20 AT OUTLET 
PARTIAL PRESSURE OF C02 AT OUTLET 
GAS TEMPERATURE AT OUTLET 
DEW POINT TEMPERATURE AT OUTLET 
GAS TEMPERATURE INSIDE HELMET 
SKIN TEMPERATURE 

LCG COOLANT TEMPERATURE AT INLET 
LCG COOLANT TEMPERATURE AT OUTLET 
LCG COOLANT FLOW RATE 
02 ADDITION/CONSUMPTION BY MAN 
02 ADDITION/CONSUMPTION BY MAN 
C02 PRODUCTION RATE BY MAN 
GAS VOLUMETRIC FLOW RATE AT SUIT INLET 
SENSIBLE HEAT GENERATED BY MAN 
RESPIRATORY LATENT HEAT GENERATED BY MAN 
LATENT HEAT GENERATED BY SWEATING 
TOTAL LATENT HEAT GENERATED BY MAN 
HEAT LOAD FROM MAN TO LCG 

HEAT LOAD FROM MAN PLUS SUIT ENVIRONMENTAL TO LCG 

LIOH BED TEMPERATURE 

GAS TEMPERATURE INSIDE LIOH CANISTER 

LIOH CANISTER WALL TEMPERATURE 

H20 COOLANT TEMPERATURE AT INLET 

H20 COOLANT TEMPERATURE AT OUTLET 

TOTAL GAS PRESSURE AT LIOH OUTLET/FAN INLET 

SENSIBLE HEAT GENERATED BY LIOH 

LATENT HEAT GENERATED BY LIOH 

LIOH SENSIBLE HEAT ABSORBED BY GAS LOOP 

LIOH SENSIBLE HEAT ABSORBED BY H20 COOLANT LOOP 

H20 FLOW RATE THROUGH SUBLIMATOR 

H20 FLOW RATE BYPASSING SUBLIMATOR 

FEEDWATER FLOW RATE 

H20 SEPARATED BY WATER SEPARATOR 

SLURPER VOLUMETRIC FLOW RATE 

CAS VOLUMETRIC FLOW RATE AT SUBLIMATOR OUTLET 

H20 COOLANT TEMPERATURE AT VEHICLE HX INLET 

H20 COOLANT TEMPERATURE AT VEHICLE HX OUTLET 

LATENT HEAT REMOVED FROM SUBLIMATOR VENT LOOP 

SENSIBLE HEAT REMOVED FROM SUBLIMATOR VENT LOOP 

TOTAL HEAT REMOVED FROM SUBLIMATOR VENT LOOP 

HEAT LOAD REMOVED FROM SUBLIMATOR H20 LOOP 

TOTAL HEAT LOAD REMOVED FROM SUBLIMATOR THROUGH SUBLIMATION 

HEAT LOAD REMOVED BY VEHICLE HX 


UNIT 

LB /HR 

LB/HR 

LB /HR 

MMHG 

MMHG 

°F 

°F 

LB/HR 

LB/HR 

LB/HR 

MMHG 

MMHG 

°F 

Op 

°p 

Op 

Op 

°F 

LB/HR 

LB/HR 

CFM 

LB/HR 

CFM 

BTU/HR 

BTU/HR 

BTU/HR 

BTU/HR 

BTU/HR 

BTU/HR 

°F 

Op 

°F 

°F 

°F 

PSIA 

BTU/HR 

BTU/HR 

BTU/HR 

BTU/HR 

LB /HR 

LB/HR 

LB/HR 

LB/HR 

CFM 

CFM 

°F 

°F 

BTU/HR 

BTU/HR 

BTU/HR 

BTU/HR 

BTU/HR 

BTU/HR 


ORIGINAL PAGE IS 
OF POOR QUALITY 


TABLE V SUITOUT PROGRAM LISTING 


10 REM "SUITOUT" -PRINTOUT FOR SUITSK2 PROGRAM WITH PLOT CAPABILI 
TV 

50 COM 01(53) ,05(18) ,03(19) ,04(8), 05 (3) , Q6( 14 ) , U1 (4) , PI (2) ,P2<16 
) ,P3(7),P4(6>,R4 

30 COM W1 ( 19 ) , W2( 1 ) , W3 ( 1 ) , T1 (54 ) , T2 ( 7 ) , T3 ( 12) , T4 (59) , T5 (5 ) , TO (5 ) 
, VI ( 17) , V2(8) , V3( 16) ,04(7) 

40 DEFFN2 ( X ) =35, 15789+24. 592588*X+2, 1 182069*Xt2-. 3414474*Xt 3+ . 15 

741642*X*4-.03i329585*X*5+.0038658282*X*6-.00024901784*X*7+.0000 

068401 559*Xt 8 

50 P3(2)=P3(1)*51.73 

GO P3(6)=P3<6>*51.73 

70 V5=V 1(5) *R5*T1 ( 4 ) /P 1 < 2 ) / 144/60 

SO 04= VI ( l)*R4*TK4)/Pl(2)/144/60 

90 Vl*=V2(8)*48.3*Tl(4)/Pl<2)/144/60 

100 T 1 ( 4 ) =T 1(4) -459 . 6 

110 T 1 =T3 ( 5 ) ~ ( ( T3 ( 5 ) ~T1 (4) ) /EXP (US/ (VI ( 5)*. 25) ) ) 

15.0 T1 ( 6) =T1 (6 ) -459.6 
130 T 1 ( 7 ) =T 1 ( 7 ) -459 » 6 
140 T1 (5)=T1 (5) -459. 6 
150 T5( 1 )=T5( 1 ) -459. 6 
160 T5(2) =T5(2) -459. 6 
170 T6( 1 ) =T6( 1 ) -459. 6 
ISO V5=» 0105*C5 
190 08=02(8) +03(8) 

200 06=01(8) -08-03(6) -03(18) 

510 P 3 ( 7 ) =P3 ( 7 ) #5 1 » 73 
220 02 ( IS) =Q2(7)-Q3( 19) 

530 R=LDG ( 10*P2 < 16 ) ) : T2 ( 6 ) =FN2 ( R ) 

540 R=L0G(1Q*P2(5) >:T2(5)=FN2(R) 

250 T2 ( 7 ) =T2 ( 7 ) -459 . 6 
260 W5= V4 ( 7 ) - V4 ( 4 ) 

570 R=L0G(10*P2(7) ) : T2<7) =FN2<R ) 

580 IF K = 1 THEN 290: OK 13) =0 
290 GOTO 580 

300 SELECT PRINT 215(132) 

310 PRINT : PR I NT 

320 PRINT "PROGRAM OUTPUT DATA: " : PR I NT 
330 PRINT "SUIT CONDITIONS: " : PRINT 

340 PR I NT " WTOTI , WH20 I „ WC02 I , PH20 I , PC02 I , TG AS I , TDP I = " ; V 1 ( 5 ) ; V4 ( 4 
) ‘ V3 ( 1 ) ' P2 ( 5 ) : P 3 ( 2 ) 8 T 1 ( 4 ) * T2 ( 5 ) 

35 0 PRINT " WTOTO, WH200, WC020, PH200, PC020, TGASO, TDPO= " ; VI ( 6 ) ; V4 ( 6 
) ; V‘3 ( 6 ) ; P2 ( 1 6 ) ;; P 3 ( 6 ) ; T*1 ( 6 ) ; T2 < 6 ) 

360 P R I NT " THEL.. , TSK I N , TLCG I , TLCGO= " ; T 1 ; T3 ( 2 ) ; T4 ( 1 6 ) ; T4 ( 6 ) 

370 PRINT " WSUIT , W02, V02, WC02, VSI =" ; W1 (2) ; V2(8) ; VI ; V3(8 ) ; V5 
380 PRINT " QMANS , QMANRL , QSW , QMANTL , QLCG , QLCGT= " ; 03 '( 6 ) ; 02 ( 8 ) ; 03 ( 8 
) ; 08: 0.6; 01 ( 18) 

390 PRINT : PRINT "LIOH CONDITIONS: " : PRINT 

. 400 PR I NT " WTOTO, WH200, WC020, PH200, PC020, TGASO, TDPO= " ; VI ( 7 ) ; V4 ( 7 


TABLE V CONT 


) : V3 < 7 > ; ps ( 7 > ; P3 < 7 > ; ti ( 7 > : ia(7) 

410 PRINT "TLIOH, TGQ, TCASE, THSQI ,THSOO, PTOTO-" ; T5( 1 ) ; T5(S) J T6< 1 ) 
:t4(S);t4(19);pki> 

420 PR I NT " QLS , QLL , QLSG , QLH2Q= " ; Q2 < 7 > ; Q2 ( 1 7 ) ; Q3 < 1 9 ) ; Q2 < 1 3 > 

430 PRINT : PRINT "SUBLIMATOR CONDITIONS: "SPRINT 

440 PR I NT " WSUBL, WB YPS , WFEED , WSEP , VSLURP , VSUBO= " ; W1 ( 8 ) ; W1 ( IS) J W1 

(io)-:w5:vs;V4 

450 PR I NT " TG AS I , THSO I , TH200 , TVHX I , TVHXO " ; T 1 ( 2 ) ; T4 < 1 3 > : T4 ( 1 4 ) ; T 


4(8 ) ; T4 ( 18) 

460 PRINT " QSUBL, QSUBS, QSUBVT, QWTOT, QTSUB , QVHXT= " ; Q1 ( 7 ) ; Q3 ( 3 ) ; Q1 
( 3 ) “ Q1 ( 23 ) ‘ Q1 ( 33 ) £ 01 (13) 

470" PRINT : PRINT "REVISED PROGRAM INPUT DATA: " : PRINT 

480 PRINT " QMETT, QENV, QFAN, QPUMP= n ; Q1 ( 8 ) ; G5 ( S ) ; Q6 ( 13 > ; Q6 ( 14 ) 

490 PRINT "PT0T„ PC02, CFM, WH20, HAST-" ; PI (2) ; P3( 1 ) ; C5; W1 ( 1 ) : U1 
500 PRINT " TWI , WCPW, CLIOH, CL2=" ; T4 ( 3 ) ; W2( 1 > JCI ;OB . 

510 PRINT "CL3, CL4, RH=" ; C3; C4; HI 

520 PR I NT " TF I , TSUBO , XKEY , UAX , WL I OH , AP , UAW= " ; T 1 < 1 3 ) ; T 1 < 5 ) ; K ; U2; W 
1(19);A1;U3 

530 PRINT "RFINn HAG AS, COP, UAHEL.= " " R 1 5 U4; P4( 1 ) ; U5 

540 PRINT “ TSK I NT, QDCMEN, QDCMEL., QEVCS, QPLSEN, I NO " ; T3( 1 ) ; Q1 < 9 ) ; Q 

kio);Q5(2);Q3(io): i 

550 SELECT PRINT 005(64) : INPUT "DO YOU WISH TO PLOT CHART (0=N0, 
1==YES> " , L. 

560 IF L=0 THEN 570:G0T0 1180 
570 END 

580 N-V1( 5 ) : GOSUB 2370: VI ( 5 > =N4 
590 N-V4 ( 4 ) : GOSUB 2370: V4 ( 4 ) =N4 
600 N~V3 ( 1 ) : GOSUB 2370: V3 ( 1 ) =N4 
610 N=P3 ( 2 ) : GOSUB 2450: P3 < 2 ) =N4 
620 N=P2(5> : GOSUB 2450: P2 ( 5 > =N4 
630 N=T 1(4): GOSUB 2330: T1 ( 4 ) =N4 
640 N-Tl; GOSUB 2330:T1=N4 
650 N==V1 (€■) : GOSUB 2370: VI (6) =N4 


660 N=V4(6) : GOSUB 2370: V4(6> =N4 
670 N=V3(6): GOSUB 2370: V3<6) =N4 
680 N=P3(6): GOSUB 2450: P3(6) =N4 


690 N=P2( 16): GOSUB 
700 N=T1 (6 ) : GOSUB 
710 N-Vl ( 7 ) : GOSUB 
720 N=V4<7>: GOSUB 
730 N-V3(7) : GOSUB 
740 N=P3(7) : GOSUB 


2450: P2( 16) =N4 
2330: T1 (6) =N4 
2370: VI ( 7 ) =N4 
2370: V4(7) =N4 
2370: V3(7) ~N4 
2450 : P 3 ( 7 > -N4 


750 N=P2(7) : GOSUB 2450: P2(7) =N4 


760 N=Tl (7) : GOSUB 2330: T1 (7) =N4 
770 N=T 1(2): GOSUB 2330: T1 (2) =N4 


780 N=T5( 1 ) : GOSUB 2330: T5( 1 ) ==N4 
790 N=T5<2) : GOSUB 2330: T5(2) ==N4 
800 N=T6C 1 ) : GOSUB 2330: T6( 1 ) =N4 
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TABLE V CONT 


-310 W-T4 ( 19 ) I G0SLJ8 0.330: T4CI.9) =N4 
850 N=T4 ( 16 ) : GOSUB 5370: T4 ( 16 ) =N4 
830 N=T4(6) : GOSUB 5370: T4 ( 6 ) =N4 
840 N=T4( 13): GOSUB 5370: T4( 13) =N4 
350 N=W1 <8> : GOSUB 5370: W1 (8 ) =N4 
360 N=T4< 18): GOSUB 5370: T4 ( 18 ) =N4 
370 1x1=14 < 14): GOSUB 5370: T4Q4 ) =N4 
380 N=W1 (5)5 GOSUB 5370: W1 (5) =N4 
890 N=Q1 (7) : GOSUB 5330: Q1 (7 ) =N4 
900 N=Q3( 3) : GOSUB 5330: G3 ( 3 ) =N4 
910 N=Q1 (3) : GOSUB 5330: 01 < 3 ) =N4 
950 M=Q 1(53): GOSUB 5330: Q1 (53) =N4 
930 N=Q1 ( 18) : GOSUB 5330: Q1 ( 18) =N4 
940 N=Q3(6): GOSUB 5330: 03(6) =N4 
950 N=Q5(7): GOSUB 5.330:05(7) =N4 
960 N=Q3( 19) : GOSUB 5330: Q3( 19) =N4 
970 N«G5( 17): GOSUB 5330: 05 ( 17 ) =N4 
980 N=G3 (8) : GOSUB 5330: 03(8) =N4 
990 N=Q5( 18) : GOSUB 5330: 05(18) =N4 
1000 N=V5: GOSUB 5370:V5=N4 
1010 N=V4: GOSUB 5370:V4=N4 
1050 N=Vl: GOSUB 5450:V1=N4 
1030 N=T5(5) : GOSUB 5330: T5 ( 5 > =N4 
1040 N=T5(6): GOSUB 5330: T5 ( 6 ) =N4 
1050 N=T5 ( 7 ) : GOSUB 5330: T5<7)=N4 
1060 N=V5(8):GDSUB 5410: V5 ( 8 ) =N4 
1070 N=W5: GOSUB 5450:W5=N4 
1080 N=T4(5) : GOSUB 5330: T4(5) =N4 
1090 N=Q8: GOSUB 5330:Q8=N4 
1100 N=Q6: GOSUB 5330:Q6=N4 
1110 N=Q1( 33): GOSUB 5330: G1 < 33 ) =N4 
1150 N=G1< 13): GOSUB 5330: Q1 ( 13) =N4 
1130 N=Wl(ia>i GOSUB 5370: W1 ( 15) =N4 
1140 N=VS: GOSUB 5370:'V5»N4 
1150 N-Wl (10): GOSUB 5410: W1 ( 10 ) =N4 
1160 N=P1(1):G0SUB 5450: P 1 ( 1 ) =N4 
1170 GOTO 300 

1180 N=P 3 ( 5 ) : GOSUB 5370: P3 ( 5 ) =N4 
1190 N=P3(6) : GOSUB 5370: P3(6) =N4 
1500 N=T4 ( 16) : GOSUB 5330: T4( 16) =N4 
1510 N=T4(6): GOSUB 5330: T4 ( 6 ) =N4 
1550 N=T4( 18): GOSUB 5330: T4( 18)=N4 
1530 N=W1 ( 8 ) : GOSUB 5330: W1 (8) =N4 
1540 N=T4( 13) : GOSUB 5330: T4( 13) =N4 
1550 N=T4( 14): GOSUB 5330: T4 ( 14 ) =N4 
1560 N=W1 (5) : GOSUB 5330: W1 (5) =N4 
1570 N=Q1 (7) : GOSUB 5590: G1 (7) =N4 
1580 N=Q3( 3) : GOSUB 5590: Q3( 3) =N4 
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1590 N=Q1 (3) ’ GOSUB 5.590: Q1 ( 3 ) =N4 
1300 N=Q1< 53): GOSUB 5590: Q1 (53) =N4 
1310 N=G1 ( IS) : GOSUB 5590: Q1 < 18 ) =N4 
1350 N=Q3<6) : GOSUB 5590: Q3(6) =N4 
1330 N=Q5(7) : GOSUB 5590: Q5(7) =N4 
1340 N=Q3( 19): GOSUB 5590: Q3< 19 ) ==N4 
1350 N=Q5( 17): GOSUB 5590: Q5< 17) =N4 
1360 N=G3(8) : GOSUB 5590: Q3< 8 ) =N4 
1 370 N=Q5 (IS): GOSUB 5590 : G5 ( 1 3 ) =N4 
1380 N=W5: GOSUB 2410:W5=N4 
1390 N=Q8: GOSUB 2590:Q8=N4 
1400 N=Q6: GOSUB 5S9Q:Q6=N4 
1410 N~Q1 ( 33) : GOSUB 5590: Q1 ( 33) =N4 
1450 N=Q1( 13): GOSUB 2290:01 ( 13) =N4 
1430 N=W1< 15): GOSUB 2330: Ml ( 12) «N4 
1440 N-P 1(1): GOSUB 2370: PI ( 1 )=N4 
1450 N=V3 ( 8 ) : GOSUB 2410: V3(8) ~N4 
1460 N=Q5 ( 8 ) : GOSUB 2290: Q2(8) -N4 

1470 PRINT : PRINT : PRINT "IS PLOTTER POWER ON?" i PRINT "IS PLOT P 
APER LOADED?": PRINT "IS CHART SWITCH ON HOLD?": PRINT "IS PEN IN 
DOWN POSITION AND CAP OFF?" 

1480 PRINT "ARE SCALE LIMITS SET (0,0 BY 10, 10 ) ? " : INPUT "IS ANSW 
ER TO ALL ABOVE QUESTIONS YES? (0=N0, 1=YES) " , L: IF LOl THEN 1470 
1490 GOTO 1 500 

1500 REM PLOT ROUTINE FOR SUITSKR 
1510 SELECT PLOT 414: SELECT PRINT 414 
1550 PLOT <5, ,C>, <24, , S> 

,iB30 PLOT <530, 775, U>: PRINT "SHUTTLE EMU FLOW CHART": PLOT <,,R> 

1 540 PLOT < 1 , , C > , < 15,, S> 

1550 PLOT <5.40, 738, U>: PRINT "CONDITIONS: MODE-":IF K = 1 THEN 1560 
: PRINT "EV":PLOT < 1 84, -50, u>: GOTO 1570 
1560 PRINT " IV" ; PLOT <-84, -50, U> 

1570 CONVERT Q1 (8) TO 21$, <-####): PRINT "METABOLIC RATE- " ," Zl$: PL 
OT <15, ,U>: PRINT "BTU/HR" : PLOT <-355,-50,U> 

1580 CONVERT Q5(8) +Q3( 10)+Q1 (9) TO 25$, (-####): PR INT "ENVIRONMEM 
T- " ; 25$: PLOT <12, , u>: PRINT "BTU/HR" : PLOT <, , R> 

1590 PLOT <400, 660, U>: CONVERT PI (5) TO Al$, ’(-##,##): PRINT "P = ";A 
1$:PL_QT <-96, -20, U> 

1600 CONVERT P3(5) TO A2$, (-#.##): PRINT " PC05=--" ; A5$: PLOT <-150,- 
20, U> 

1610 CONVERT T 1 ( 4 ) TO A3$, <-##.#): PRINT "T=" ; A3$: PLOT <-84,-20,U 

1650 CONVERT T5(5) TO A4$, (-##.#>: PRINT " DP - " ; A4$ : PLOT < -96 , -20 , 
U> 

1630 CONVERT VS TO A5$, (-#.##): PR INT " V= " ? A5$: PLOT <-84,-50,U> 
1640 CONVERT VI (S') TO A6$, (-##,##)■ PRINT " MTOT - " ; A6$ : PLOT <,,R> 
1650 PLOT <400, 460, U>: CONVERT P3(6) TO Bl$, (-#.##): PRINT "PC05=" 
J B1$S PLOT <-150, -20, U> 
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1 660 CONVERT T 1(6) TO B2* , ( -## „ # ) : PR I NT " T= " ; BE* : PLOT <-84 ,-£0,0 

1670 CONVERT T3(6) TO B3*, (-##„#) I PRINT "DP-" ; B3*: PLOT <, ,R> 

1660 PLOT <520, 440, U>: CONVERT 03(16) TO Cl*, (-###> : PR INT " QSEN=" 
; C .1 $: PLOT <-108, -£0, u> 

1690 CONVERT 03(17) TO C2*, (-###) SPRINT 11 QLAT=" S CE*S PLOT <,,R> 
1700 PLOT <570, 385, U>: CONVERT T 1 ( 7 ) TO D1 *,<-###„#): PRINT ;, T-'^D 
1$SPL0T <-96, -EQ, U> 

1710 CONVERT T£(7) TO DE$, (-##.#) S PRINT “ DP=" ; D£*S PLOT <,,R> 

1720 PLOT <740, 440, U>: CONVERT C5 TO El $,(-#.##) I PRINT "V= : ‘ ; E1*;P 
LOT <-84, -20, U> 

1730 CONVERT Pl(l) TO EE*, <-##„##) S PRINT "P=" ; EE*S PLOT <,,R> 

1740 PLOT <820, 370, U>: CONVERT 06(13) TO FI*, <-##) SPRINT "Q=":;F1* 
: PLOT < , , R > 

17S0 PLOT <857, 440, U>! CONVERT Tl<2> TO Gl$, (-###.#): PRINT "T-";G 
1*S PLOT <-96,-20,U> 

1760 CONVERT T2(7) TO GE*, (-##„#): PRINT " DP- " ; GE$: PLOT <,,R> 

1770 PLOT <983, 500, U> SPRINT "VENTS": PLOT <-60,-£0,U> 

1780 CONVERT Q1 ( 3 ) TO Hl$, (-###) 5 PRINT " QTOT= H1*S PLOT <-108, -2 
Q,U> 

1790 CONVERT Q1 (7) TO HE* , ( -### ) S PRINT " QLAT== " ; HE*: PLOT <-108, -2 
0,U> 

1 800 CONVERT Q3 ( 3 ) TO H3$ , ( -### ) : PR I NT " Q8EN- " S H3* S PLOT < , , R> 
1810 PLOT <790, 530, U>S CONVERT VE TO 11$, <-#.##> SPRINT ,, V* i, ;Il$:P 
LOT <, ,R> 

1820 PLOT <790, 660, U>S CONVERT Wl(10) TO Jl$, (-#.##) SPRINT "MH20= 
" Ml PLOT <,,R> 

1830 PLOT <710, 590, U>S CONVERT W5 TO Kl$, <-#.###) SPRINT "MH£0=" ;K 
1$: PLOT <„,R> 

1840 PLOT <880, 760, U>S CONVERT T 1(4) TO Ll$, (-##.#) SPRINT "T=";L1 
$ S PLOT < -84 , -30, LJ> 

18B0 PRINT " DP=" ; LI*: PLOT <-96,-20,U> 

1860 CONVERT V4 TO L3$, (-#.##) s PRINT " V=" ; LE*: PLOT <-S4,-20,U> 
1870 CONVERT Vl(5)-V3(8) TO L3$, <-##.##) S PRINT " M=" ; L3$S PLOT <, , 
R> 

1880 PLOT <540, 615, U>S CONVERT VE(8) TO Ml$, (-#.###) SPRINT "MOE=" 
SPLOT <-48, -£0,U>S PRINT Ml*: PLOT <,,R> 

1890 PLOT <300, 380, U>S CONVERT T4(6) TO N1 $,(-###.#) SPRINT "T=“ ;N 
1$SPL0T <,,R> 

1900 PLOT <470 , 23 l,Li>: CONVERT 03(10) TO 01$, (-###): PRINT "Q=";01 
* SPLOT <t72, -25, U> 

1910 CONVERT 05(E) TO 02$, <-###): PRINT "Q=" ; 02*: PLOT <,,R> 

19E0 PLOT <620, 196, U> 5 CONVERT 06(14) TO Pl$, (-###) SPRINT "Q=";P1 
*: PLOT < , , R > 

1930 PLOT <600, 330, U>S CONVERT T4(£) TO 01*, (-###.#): PRINT "T= ,, S Q 
1$ SPLOT <-192,,U> 

1940 CONVERT 02(18) TO Rl$, (-###) S PRINT "Q=" ; Rl$s PLOT <-19£,,U> 
1950 CONVERT WK19) TO SI*, (-##): PRINT " M-'" ; SI*: PLOT <-60,-20,U> 
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.'I960 

1970 


CONVERT T4 ( 19 ) 
PLOT <710,250, 


TO S2$, <-###. ■#>: PRINT " T~ " S2$; PLOT <,,«;■ 
U>: CONVERT T4CI.B) TO Tl$, (-###..#): PR I NT " T - 


.1 


Tl$; PLOT < , , R > 

1980 PLOT <730, 330, U>: CONVERT 01(10) TO 01$, (-###): PRINT " QELEC= 
" ; Ul$: PLOT <-120, -20, u> 


1990 CONVERT Q1 (9) TO 02$, (-###): PRINT "QENVIR=" ; 02$ I PLOT <, ,R> 
2000 PLOT <840,250, 0>: CONVERT T4(18> TO 03$, <-###.#): PRINT "T=" 
03$; PLOT < , , R > 


2010 PLOT <984, 330, 0>: IF K-l THEN 2020; CONVERT 01(33) TO 04$, ( -# 

###): PRINT " QTC.lT- 11 ; 04$: PLOT <,,R>:GOTO 2030 

2020 CONVERT V TO 04$, <•-#); PRINT " QTOT-" 5 04$: PLOT <,,R> 

2030 PLOT <988,200, U> : CONVERT Q1 ( 23 ) TO 05$, ( -#### ) : PRINT " QH20= 
" : PLOT <-60, -20, 0>: PRINT 05$; PLOT <-60, -20, 0> 

2040 CONVERT Wl(8) TO 06$, (-###.#): PRINT “M*" J 06$; PLOT <-96,-20, 


2050 CONVERT T4C14) TO 07$, (-##,#) C PR INI "T~" ; 07$; PLOT <,,R> 

2060 PLOT <830, 140, 0>; IF K = i THEN 2070; CONVERT V TO Vl$, (.-#.#): P 
PINT " " ; V 1 $ : PLOT <,,R>:GOTO 2080 

2070 CONVERT W 1 ( 8 ) TO Vl$, (-###.#): PR I NT " M= " : Vl$; PLOT <,,R> 

2080 PLOT <630, 170, 0>; CONVERT W1CL2) TO V2$, (-###.#): PRINT "M= M ; 
V2$; PLOT <,,R> 

2090 PLOT <898, 75, U>: IF K-l THEN 2100: CONVERT V TO V3$, <-#.#): PR 
I NT " T=" ; V3$; PLOT <-72, -25, u>: GOTO 2.110 

2100 CONVERT T4 (8 ) TO V3$„ ( -##. # ) : PRINT " T - " ; V3$; PLOT <-84, -25, U 
> ; GOTO 2.120 

2110 CONVERT V TO V4$, ( -# ) : PR I NT " MCP = " ; V4$; PLOT <,,R>:GOTO 2130 
2120 CONVERT W2(l) TO V4$, (-###); PRINT " MCP - " ; V4$ : PLOT <,,R> 

2130 PLOT <780, 80, U> S' IF K==l THEN 2140: CONVERT "V TO V5$, ( -#) : PRIN 
T "G=" : V5$; PLOT <, , R>: GOTO 2150 

2140 CONVERT Q1 (13) TO V5$, ( -#### ) : PR I NT " Q= " ; V5$; PLOT < , , R> 

2150 PLOT <280, 120, U>: CONVERT W 1 ( 2 ) TO V6$, (-###.#): PRINT "M=";V 
6$; PLOT <-96, -30, U> 

2160 CONVERT T4 ( 16 ) TO V7$, ( -###,. # ) : PRINT " T-" ”, V7$; PLOT < , , R> 
2170 PLOT <110, 605, U>; CONVERT T1 TO Wl$, (-##.#) S PRINT "T=";PLOT 
<-18, , U> : PR INT' wi$; PLOT < , , R> 

2180 PLOT <220, 660, U>: PRINT "RESP; “ : PLOT <-60,-20,U> 

2190 CONVERT 02(8) TO W2$, (-###): PRINT " QLAT - " ; W2$ ; PLOT <-108,-2 

o,u> 

2200 CONVERT V3<8) TO W3$, <-#.###): PRINT " MC02- " ; W3$ S PLOT <,,R> 
2210 PLOT <94, 510, U>* PR INT "VENT-" : PLOT <-88,-20,U> 

3220 PRINT "SUIT TOT! "S PLOT <-108, -20, LJ> 

2230 CONVERT 08 TO W4$» (-###); PRINT " QLAT ~ " ; W4$; PLOT <-108, -20,11 

2240 CONVERT 03(6) TO W5$, (-###): PRINT "O.SEN-" ; W5$: PLOT <,,R> 
2250 PLOT <90, 50, U>S PRINT " LOG CONVERT 06 TO W6$, (-####): PRINT 
"QMET~" ; W6$ ; PLOT <-.120, -20, U> 

2260 CONVERT 05(8) TO W7$, (-####); PR INT " QENVIR= " ; W7$: PLOT <-144 
, -20, LJ> 
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5870 CONVERT QIUSTtS W8$, (-####): PRINT " QTOT= " ; WS$ ; PLOT <,,R 
2280 STOP 
2.290 N3=INTCN) 

2300 D=N~N3MF D>=„5 THEN 2310; GOTO 2320 

2310 N4=N3-f-lCREETURN 

2320 N4=N3; RETURN 

2330 N1=N*10:N2-XNT<N1 ) :N3=N5/10 

2340 D=N-N3: IF D>=„05 THEN 2350: GOTO 2300 

2350 N4-N3+ „ 1 : RETURN 

2360 N4»N3: RETURN 

2370 N1=N*100;N2=INT(N1) :N3=N2/100 
2380 D=N-N3“.IF D>*.0Q5 THEN 2300: GOTO 2400 
2390 N4=N3+» 01 : RETURN 
2400 N4=N3: RETURN 

2410 N1=N*1000* N2-INT (N1 ) :N3=N5/1000 
2420 D=N-N3: IF D>=.0005 THEN 2430; GOTO 2440 
2430 N4=N3+ . 001 : RETURN 
2440 N4«N3: RETURN 

2450 N1=N*1QQ0G:N2"INT<N1 ) :N3=N2/10000 
2460 D-N”N3: IF D>*. 00005 THEN 2470: GOTO 2480 
2470 N4=N3+ „ 000 1 : RETURN 
2480 N4=N3: RETURN 
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LIOHPPC COMPUTER PROGRAM 


File Name "LIOHPPC" 

Abstract "LIOHPPC" plots the partial pressure of C02 at the upper 
and lower LIOH canister outlets and at the ARS inlet. 

The program also provides an optional printout of the 
plotted data. Figures 5 and 0 provide examples of the 
plot and its corresponding data printout. The program 
is designed for use with the Wang 2200 series computer 
system. 

Program Description 

The program supplies detail instructions on the CRT while it is being 
run. Initially, the user chooses between a plot of the data and a 
printout of the data. If plots are required, the user is asked whether 
the data will be keyed in or if it has been previously stored on tape. 
If it has been previously stored on tapes, the program will instruct 
the user to load the appropriate tape and file number and the data 
will be loaded. If the data is to be keyed in, the CRT will ask for 
the dat.e of the test, the type of case (3.75 man case, 7.25 man case, 
etc.), and the value of x at the origin. The program then requests 
the test data points. The user inputs the time in minutes, the par- 
tial pressure of CO 2 at the ARS inlet in mmHg, the partial pressure 
of CO 2 at the upper canister outlet in mmHg, and the partial pressure 
of CO 2 at the lower canister outlet in mmHg. A total of 150 sets of 
data points can be keyed in by the user. Before the data is stored 
onto tape, the program allows corrections to be made to the data. 

This can be accomplished by inserting line numbers 261 through 269 
into the program with revised data and executing these line numbers. 

The program then instructs the user to place a data tape into the 
tape drive 10B, asks which file number is to be used, and saves the 
data. 

Next, the program asks whether plots are to be made. If the user 
answers no, the program is ended. But if plots are desired, the CRT 
asks whether a 900 minute or 1800 minute time range is desired. This 
option allows the user to choose between a 900 or an 1800 minute time 
range along the x axis. This was necessary because some of the cases 
had longer times between canister changeouts than other cases. When 
all of the required plot parameters are defined, the program utilizes 
the generalized plot routine as a subroutine to produce the plots. 

Next the program instructs the user to load the graph paper onto the 
flatbed plotter and then produces the graph exemplified by the 7 man 
case shown in Figure 5. The label on the plot is optional and if 
desired, the user must position the label where it will not interfere 
with the graph. 
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> 

If initially a printout of the results was desired instead of plots, 

the user is instructed to load the appropriate data tape onto the tape 

drive unit 10A and specify which file is to be printed out. The pro- | 

gram then supplies the print out as exemplified by the 7 man case shown f 

in Figure 6. A listing of the program "LIOHPPC" is provided in Table 

VI. The data used to generate the plots has been saved and retained 1 

on tapes for possible further analysis. A discussion of the RSECS I 

LiOH test results based upon the plots provided by this program can be 

found in Memorandum EC2-77-185 dated 11-8-77. . 






i 


.L . 


ORIGINAL PAGE IS 
OF POOR QUALITY 







37 


FIGURE 6 


LIOHPPC PROGRAM EXAMPLE DATA TABULATION 


R S I • C S ( •' P C02 D A T A ! F I L. E N U M B E R 


1 


T E >3 T D i-y T E -- 4 / 2 0 / 7 7 

TYPE op CASE --7 MAN CASE 

VALUE OF X AT THE ORIGIN-- 3300 

TUlAL NUMBER OF DATA POINTS™ 108 


I 

1 

A.. 

3 

A 

5 

6 

7 

8 
9 

TO 

II 
12 

13 

14 
TS 
16 

17 

18 

19 

20 


24 

26 

26 

27 

28 

29 

30 

31 

32 

33 

34 
36 

36 

37 

38 

39 


TIME 
Q ( I ) 

3329 

3330 

3331 

3332 
3336 
3340 
3346 
3360 
3366 
3360 
3366 
3370 
3376 
3380 
3390 
3400 
3411 
3420 
3428 
3464 
3460 
3468 
3482 
3490 
3495 
3518 
3527 
3630 
3540 
3550 
3655 
3562 
3570 
3580 
3590 
3600 
3610 
3620 
3630 


ARSIN 
H ( I > 
7. 57 
7 1 , 38 

7.07 
6. 79 
5., 92 
4 u 8 4 
4u 1 

3 . 61 
3.2 
2 . 87 

2.7 
2,47 

39 
27 
1 6 
09 
1 

1 5 
17 

•U li iUi w 
2.28 
2 „ 34 
2.49 
2 . 6 
2.66 
2 . 88 
3.01 
3,04 
3. 19 
3.24 
3. 33 
3, 37 
3.6 
3. 58 
3.67 
3.83 

3. 93 

4. 07 

4 . 1 6 


UP CAN 
P < I ) 
6.92 
6, 87 
.21 


LOWCAN 
Q 1(1) 
6.98 
, 33 
.39 


. 17 

,37 

.28 

.68 

» 46 

, 75 

» 5 

.79 

.65 

. 79 

„ 69 

. 78 

, 59 

.84 

, 6 

, 9 

,66 

-88 

» 66 

-88 

,73 

.91 

,78 

.96 

, 86 

1 . 03 

, 96 

1.06 

1 . 03 

1.18 

1 . 1 2 

1 . 24 

1, 43 

1.48 

1.49 

1.48 

1 . 57 

1 . 54 

1 . 7 

1.69 

1 „ 78 

1. 72 

1.85 

1.85 

2 . 09 

2.06 

2.24 

9 9 


':> t: 

« 4 .W 

2 . 19 

2.41 

2.29 

2 . 53 

2.44 

■"» C7 

2.49 

2 , 6 4 

2 5 ^ 

uQ 

7 7 

aL n / O 

2 . 61 

2.88 

2 . 69 

2 . 97 

2 . 86 

3 , 11 

2.99 

3.21 

3 . 08 

3 , 3 

3.23 

3 . A 1 

3.35 
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40 

41 

42 
4 3 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 
SB 

59 

60 
61 
62 
63 
6 4 
6 S 
66 
A / 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 


3640 

4 . 29 

3650 

4. 38 

3660 

4. 52 

3670 

An 64 

3680 

4.67 

3690 

4.78 

3700 

4* 93 

3710 

5-01 

3720 

c -i c 

W a 4. 

3730 

5-27 

3740 

5.36 

3 7 SO 

5 - 55 

3759 

S - 62 

3770 

5. 76 

3780 

5-89 

3790 

6-0 4 

3800 

6- 18 

3810 

6-28 

3819 

6 - 39 

3830 

6- 56 

3840 

6„ 7 

3848 

6-84 

3860 

6.97 

3870 

7. 1 

3877 

7. 18 

3884 

7- 32 

3890 

7-4 

3900 

7-61 

3909 

7. 71 

3915 

7. 75 

3920 

7.89 

3930 

8-03 

3940 

8-21 

3950 

8-43 

3960 

8.55 

3970 

8-67 

3980 

8.84 

3990 

9,04 

3998 

9. 17 

4009 

9-28 

4013 

9-36 

4015 

9. 06 

4016 

8-84 

4018 

8.55 

4020 

8.27 

4023 

7-87 

4027 

7. 5 


3-54 

3.41 

3 » 65 

3.53 

3-75 

3-66 

3- 9 

3.72 

3-94 

3.81 

4. 07 

3.93 

4- 19 

4-08 

4-26 

4- 17 

4-42 

4-26 

4- 55 

4. 42 

4. 58 

4- 53 

4. 7 

4. 59 

4-83 

4.67 

5. 1 

4. 83 

5-23 

4.95 

5.43 

5. 16 

5-57. 

5.24 

5.67 

5. 41 

5.87 

5.52 

6- 07 

5.74 

6. 19 

5.87 

6- 35 

5.99 

6- 8 

6. 17 

6 - 6 5 

6.24 

6-82 

6.49 

6.89 

6-65 

7.06 

6.66 

7.21 

6.81 

7-36 

6-94 

7.44 

7.08 

7- 58 

7- 13 

7- 68 

7-29 

7-83 

7.53 

8- 1 

7.62 

8.23 

7.77 

8.37 

7.98 

8- 53 

8- IS 

8.84 

8.33 

8.86 

8.47 

9. 11 

8. 64 

9.38 

3. 72 

8. 9 

2. 94 

8. 7 

3.04 

8. 4 

2.83 

8- 05 

3. 05 

7-67 

3. 12 

7- 27 

3. 18 


E 



87 

4030 

7 

| 

88 

4035 

6 

f 

89 

4040 

6 

1 

90 

4050 

5 


91 

4060 

5 


92 

4070 

5 


93 

4080 

4 


94 

4086 

4 


95 

4093 

4 


96 

4100 

4 

-#■ ■ 

97 

4110 

4 


98 

4119 

4 


99 

4130 

4 


100 

4140 

4 


101 

4148 

5 

1 

\. 

102 

4155 


; 

103 

4160 

5 


104 

4170 

5 


105 

4180 

5 

1 

106 

4197 

5 

( : 

107 

4200 

5 

j 

! 

108 

4210 

5 


5 . 83 


6.92 

6 . 5 5 

6. 19 

S. 64 

5. 28 

4.99 

4.9 

4. 86 

4.84 

4 .77 

4.81 

4.79 

4.83 

4.95 

4.98 

5.06 

5.1 

5.23 


5. 58 
5. 74 
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TABLE VI LIOHPPC PROGRAM LISTING 



20 COM X9 ( 1 50 ) „ Y9 ( 1 SO ) , C < 1 0 > , X$25 , Y$25, P$ 1 , C$30 y D$20 
30 COM £$20, F$20, Zl, G$, H$, X$ y J$ 

40 COM GK150),H(i50),P<lS0),D(10),Qia50) 

SO E$~"LIOH TEST RESULTS" 

60 G$“" ARS INLET PC02" 

70 H$ s = " UP CAN OUT PC 02" 

80 1$-"-" LOW CAN OUT PC02" 

90 PRINT "INPUT TYPE OF RUN DESIRED" ! PRINT " 1=L0AD DATA FROM TA 
PE AND OBTAIN PRINTOUT OF DESIRED FILES" 

100 PRINT " 2*L0AD DATA (TAPE OR KEY IN) AND OBTAIN PLOTS" 

110 INPUT Q; IF Q-2 THEN 120: GOTO 1960 

120 INPUT "IS DATA TO BE LOADED FROM TAPE (0=NQ , 1=YES>",G: 

IF G*0 THEN ISO 

3.30 PRINT "LOAD DATA TAPE ON UNIT B "! INPUT "FILE NUMBER TO BE L 
OADED" , N1 ! SELECT TAPE 10B I -P9*0 : REWIND : FOR K“1 TO 50:P9=P9+1: 
NEXT K! IF N 1 * 1 THEN 140 : N N 1 ~ 1 1 SKIP NF 
140 DATA LOAD "LTEST" I DATA LOAD D$, J$, Zl, D(8) , QO , H< ) , P( ) , Q1 ( > I 

IF END THEN 300: GOTO 300 
ISO INPUT "DATE OF TEST < MTH/DAY/YR > " » D$ 

160 INPUT "TYPE OF CASE".J$ 

170 INPUT " VALUE OF X AT THE ORIGIN"', Z1 

180 PRINT : PRINT "INPUT TEST DATA (STOP BY ENTERING NEGATIVE VAL 
LIES)" : PRINT 
3.90 1*1 

200 PRINT "DATA POINT TIME, ARS y UP CAN, LOW CAN 

210 PRINT " " J I , ; : INPUT Q ( I ) , H ( I ) , P ( I ) y 01 ( I ) 

220 IF 0 ( I ) >-,0000001 THEN 230 ; I - 1 -• 1 : D ( 8 ) - 1 : GOTO 2S0 
230 1*1+1 
240 GOTO 210 

2 SO PRINT "ARE THERE ANY DATA CORRECTIONS TO BE MADE?" 

260 STOP 

270 PRINT "LOAD DATA TAPE ON UNIT 10B" i PRINT 
280 INPUT "FILE NUMBER WHERE DATA IS TO SAVED", Nil 

SELECT TAPE 1 OB: REWIND JP9*0:FOR K*1 TO SO : P9*P9+1 : NEXT K! 

IF Nl=l THEN 290: N-Ni-l : SKIP NF 

290 SELECT TAPE 10B: DATA SAVE OPEN "LTEST": DATA SAVE D$,J$ f Zl,D( 
8 ) y Q ( ) y H ( )„P( ) , Q 1 C ) : DATA SAVE END : SELECT TAPE 10A 
300 SELECT TAPE 10A: INPUT "DO YOU WISH TO PLOT DATA (0=N0,1^YES) 
" yQ9:lF Q9-0 THEN 500 

310 PRINT "DO YOU WISH TO HAVE A 900 OR 1800 MINUTE TIME RANGE?" 


320 PRINT " 1*900 MINUTE TIME RANGE" I PRINT " 2*1800 MINUTE TIME 
RANGE" 

330 INPUT Q: IF G-2 THEN 3S0 


340 K6*D(8):X0* 
0: T2“9 ! X$™" TIME 


= 100! YO*it Xl*l: Y1*1'S C1*Z1 
: IN MINUTES" : Y$=="PARTIAL 


C2*o: si=zi: S2=z 1+900: Ti* 
PRESSURE OF C02" : Ll*2 : L2 
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=1! GOSUB '00: GOTO 360 

3 SO K6«D(8> : X 0-2 00 5 Y0*2‘. Xl=l : Yl = l ; Cl-Zl : C2=0: S3. =23. 1 S2=Z:L ♦1800! T1 

-0: T2-12: X$="TIME IN MINUTES" ! Y$-" PARTIAL PRESSURE CM- - C02" !L1«2I 

L.2-1: GOSUB '00 

360 FOR R=1 TO 3 

370 IF R*3 THEN 4 SO 

380 IF R-2 THEN 420 

390 FOR I = 3. TO K6 : X9 (I ) =Q (I > I f 9 ( I ) =H (I ) : NEXT I 
400 P*=HEX ( 00 ) ; 0= 1 : K= 3. : K5=K6 5 L8=l I L9=K6 : GOSUB ' 03 
410 GOTO 470 

420 FOR 1 = 1 TO K6 1X9(1) =Q (1)1 Y9 < T. ) =P (I ) ! NEXT I 
430 P$=HEX < 01 > 1 0 = 3. 1 K=1 1 K5=K6 1 L8=l 1 L9=K6 1 GOSUB ' 03 
440 GOTO 470 

4 SO FOR I = 1 TO K6 1 X9 ( I ) =Q ( I ) 1 Y9 ( I ) =Q1 (1)1 NEXT I 
460 P$=HEX ( 02 > 1 0=1 1 K=1 1 KS=K6 1 L8=l 1 L9=K6 1 GOSUB ' 03 
470 NEXT R 

480 INPUT "00 1011 WISH TO LABEL GRAPH? (0=N0, 1=YES ) “ , P 1 IF P=0 T 
HEN 500 

490 GOSUB '04! GOSUB 1890 
500 STOP 

S10 OEFFN' 001 Q=0l GOSUB 520; PLOT <,,R>; RETURN 
S20 SELECT PRINT OOS; PRINT HEX ( 03 >1 PRINT : PRINT : 

PRINT "*** GENERATING PLOTS OF " ; E$; " ***" 

530 PRINT 1 PRINT 1 PRINT " IS PLOTTER POWER ON : PRINT " IS PLOT 
PAPER LOAOEO ? " I PRINT " IS CHART SWITCH ON HOLD ? " f PRINT " IS 
PEN IN DOWN POSITION AND CAP OFF ?" 

S40 PRINT " ARE SCALE LIMITS SET (0,0 BY 10,10) ?" 

SSO INPUT " IS ANSWER TO ALL ABODE QUESTIONS YES ? ( YES=1, N0=0) " 
, Q : IF QOl THEN 520 

S60 SELECT PLOT 414; PLOT < X , , C> f < 12., S> 5 0=0 
S 7 0 Fl = 100/X0: F2--100/Y0: GOSUB 800; PLOT <1, , C>, <, , S> 

580 RETURN 

SSO OEFFN' 03: SELECT PRINT OOS! PRINT HEX (03)! 

PLOT <„ , R>, <10Q«Xi, 100*Y1, U>, <1, . C> 

600 X4 , Y4 f E„ E3, E4 » E6, E8 , E7=0 1 X$=" ": U1 = 1 
610 K=l!X4,Y4=0 
620 K1=K: W=1 

630 IJ2=0: X=X9(K) : IF X>S2 THEN 680: Y=Y9(K) ; IF Y>9 THEN 680: GOSUB 

' 3. 25 : I F U2= 3. THEN 700 : X=X9 ( K > -C 1 ! Y= Y9 < K ) -C2 : X4, Y4=0 

640 IF K>K:L THEN 660 1 PLOT <, , R>„ <100*X1, 100*Y1 , U>, <F1*X, F2*Y, U> 

6 SO PLOT <, . D>, <„ , P$>: GOSUB 1720: GOTO 670 
660 GOSUB ' 122(X, Y, X4, Y4 ) 

670 PLOT <-X*Fi , -Y*F2, U> ! GOTO 700 
680 IF K/5=W THEN 6901 GOTO 700 
690 W=W+1 

700 K=K+ 1 : IF K>=KS+ 1 THEN 710; GOTO 630 

710 X4, Y4,E f E3,E4, E6,E8,E7=0l Ul=2: PLOT <, , R>„ <100*X1 , 100*Y1 , .US- 
720 FOR I=L8 TO L9 
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730 X«X9 < I ) -C 1 ; V "Y9 < I ) -C2 

740 IF X>?00 THEN 77 0: IF Y>9 THEN 770 

750 IF 101.8 THEN 760 : Ul = l : GOSUB f 122 ( X, Y, X4 y Y4 ) ; Ul=2 ! GOTO 780 
760 GOSUB ' 122 ( X , Y , X4 , Y4 ) ! GOTO 780 
770 U3--1 
780 NEXT I 


790 
800 
( T2' 
810 
U>, 
820 
8,30 
840 

8 SO 
860 
870 
880 
890 
900 
910 
920 
930 
940 

9 SO 
9 6 0 


PLOT <, , U> I PLOT <, , R> ; RETURN 

A 1 ~F 1 #ABS ( S 1 -C 1 ) ; A2«F 1 #ABS ( 82-C 1 ) ! B1 ~F2#ABS < T 1 ~C2 ) ; B2=F2#ABS 
-02 ) 

PLOT <, , R> v <100# XI , 100*Y1 , IJ>, <-Al, O r U>, <A1+A2, 0, D>, <~A2, ~B1 , 
<0,B1+B2,D> 

MS” ( ABB < S 1-0 1 ) + ABS < 82-01 ) ) /XO : NS- < ABS < T 1 -02 ) + ABS < T2-C2 ) ) /Y 0 
K-0 

83-81 -XO 

PLOT <-3 , - ( B 1 +B2 > , U> 

NS-INT (NS+« S ) 

FOR 13=1 TO N5 + 1 
PLOT <6,0,D>,<, ,U> 

IF I3-N5+1 THEN 900; PLOT <~6 y 100, LO- 
NE XT 13 

PLOT <•- < A 1 + 3 ) , - < B2+6 ) , U> 

HS=INT(HS+n S) 

FOR 14=1 TO MS+1 
PLOT <0, 12, D>, <, , U> 

IF I4-MS+1 THEN 960; PLOT <100,-12,U> 

NEXT 14 


970 IF LI =2 THEN 980: PLOT <-< A1+A2+24 ) , 20, U>: GOTO 990 
980 PLOT <™ < A1+A2+24 > , -32, U> 

990 FOR 1=1 TO M5 + 1 
1000 IF I>M5+1 THEN 1110 
1010 S3=S3+X0 

1020 IF ABS ( S3 ) >= 1 000 „ THEN 1040; IF ABS < S3 >>=100, THEN 1050: IF ABS 

< S3) >--10. THEN 1060; IF ABS(S3)>=1., THEN 1070 

1030 CONVERT S3 TO S3*, <-.###) 5 GOTO 1080 

1040 CONVERT S3 TO S3*,- <-####): GOTO 1080 

1050 CONVERT S3 TO 83$, (-###. ) ; GOTO 1080 

1060 CONVERT S3 TO S3*, <-##.#>! GOTO 1030 

1070 CONVERT S3 TO .S3*, <-#„#«); GOTO 1080 

1080 IF KOO THEN 1090; PLOT S3$>! GOTO 1110 

1090 IF I - 1 <> IN T ( A 1 »S)/ 1 0 0 THEN 11 00: PLOT <100, 0, U>: GOTO 1110 

1 1 00 PLOT <1 00-60 . 0 , U>, < , , S3*> 

mo k=k+i ; next ;i: 

1120 IF L2=l THEN 1130; PLOT <- < A2-+-18 ) , 0, U> ; GOTO 1140 
1130 PLOT <- ( A2+90 ) , 0, U> 

1140 IF L 1 = 1 THEN 11 SO; PLOT < 0 , ••• ( B 1 -26 ) , U> I GOTO 1160 
1150 PLOT <0, -<Bl+26> , U> 

1160 K-0 
1170 T3=T1-Y0 
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1180 FOR I2=:= 1 to N5+1 
1190 IF I2>NS+1 THEN 1300 
1200 T3=T3+Y0 

1210 IF ABS ( T3 ) >=1000. THEN 1230! IF ABS (T3 ) >=100- THEN 12401 IF ABS 

<T3)>=10« THEN 1250; IF ABS ( T3 ) >=1 » THEN 1260 

1220 CONVERT T3 TO T3$, <-,##>! GOTO 1270 

1230 CONVERT T3 TO T3*, <-##«#): GOTO 1270 

1240 CONVERT T3 TO T3*. <-«#».) 1 GOTO 1270 

1250 CONVERT T3 TO T3$, #> ! GOTO 1270 

1260 CONVERT T3 TO T3$» <■-#.#) i GOTO 1270 

1270 IF KOO THEN 1280: PLOT <, , T3$> ! GOTO 1300 

1280 IF I2-iOINT<Bl+. 5)/100 THEN 1290! PLOT <0, 100, U>! GOTO 1300 
1 2 90 P L 0 T <--48,100, U > , < , , T3 $ > 

1300 K-K+l : NEXT 12 

1310 IF X$==" " THEN 1370: PLOT <, , R>, <100#X1-A1 , 100*Y1 , U> 

1320 IF L 1 ~2 THEN 1330: PLOT <0, SO, U>; GOTO 1340 
1330 PLOT <0,-60, U> 

1340 IF A 200 THEN 1350: PLOT <A:L/5, 0, U> ! GOTO 1360 
1350 PLOT <A1+A2/S, 0, U> 

1360 PLOT <..X*> 

1370 IF Y$=" " THEN 1440 

1380 PLOT <,,R>,<100*X1, 100#Y1-B1* U> 

1390 IF L2-1 THEN 1400: PLOT <90, 0, U> I GOTO 1410 
1400 PLOT <-90 , 0 , U> 

1410 IF B2O0 THEN 1420: PLOT <0, 2*B'l/3, U> ! GOTO 1430 
1420 PLOT <0,B1+B2*2/3,U> 

1430 PLOT <0, -20, S>, <, , Y*>, <13, , S>, <, , R> 

1440 RETURN 

1450 DEFFN y 122 (U, V, X4, Y4 ) 

1460 X3-U: Y3-V 

1470 D1==X3-X4: D2=Y3 -Y4: X4-X4+D1 ; Y4-Y4+D2 
1480 IF Ul==2 THEN 1540 

1490 IF U 2 =:: 1 THEN 1520: PLOT <F1*D1 , F2*D2, U >, .<, , D> 

1500 PLOT <,,[»; IF K/S==W THEN 1510: IF K=K5 THEN 1510: PLOT <,,U> 

: GOTO 1530 

1510 GOSUB 1720: W=W+l! PLOT <, , U> : GOTO 1530 
1520 PLOT <F 1 * D 1 , U 3 , U > , < » , U > * < , , >,<, ,U>: 

E3-INT ( U3 ) : E7-INT < F1*D1 ) : GOTO 1670 
1530 E3-INT ( F2#D2 ) i E7=INT ( F1*D1 ) ! RETURN 
1540 E1=F2#Y3-INT < F2*Y3 ) 

1550 E==E 1 +E4 

1 560 P8~F2*D2+E ! P9= I NT < P8 ) 

1 570 E5-F 1 «X3- 1 NT < F 1 #X3 > 

1580 E6-ES+E8 

1590 S8=:F1#D1+E6! S9-INT<S8> 

1600 IF U2=0 THEN 1610: PLOT <, , U>, <S9, 0, U> : GOTO 1630 
1610 IF U3-1 THEN 1620: PLOT <, , D>, <S9, P9, D> : GOTO 1630 
1620 PLOT <, , U>, <89, P9, U> ! U3==0 
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,1,630 IF U2-1 THEN 1640! E3==E3+P9 
1640 E4-F2' X Y3~E3 
1650 E7-S9+E7 
1660 E8-F1' X ‘X3~E7 
1670 RETURN 

1680 DEFFN' 125; IF < Y - C 2 ) / Y 0 Y 1 < 9 „ 9 9 9 T HEN 1690! U2-1 
1690 IF <Y-C2)/Y0+Y1>,001THEN 1700 ! U2*i 
1700 IF (X-CD/X0+XK9.999THEN 1710! U2*l 
1710 RETURN 

1720 ON < UAL < PI ) + 1 > GO SUB 1730, 1740, 1750 ! RETURN 

1730 PLOT <0, 8, U>, <-7, -12, D>, <14, 0, D>, <~7, 12, [», <0, -8, U>! RETURN 

1740 PLOT <0, 6, U>, <-4, -12, D>, <10, 8, D>, <~12, 0, t», <10, -8, D>, <~4, 12 
, D>, <0,-6, U>! RETURN 

1750 PLOT <-4, 7, D>, <8, 0, D>, <-8, -14, D>, <8, 0, D>, <-4,7, D>, <-7, 4, D>, 
<0, -8, D>» <14, 8, D>, <0, -8, D>, <-7, 4, D>, <0, 0, U> ! RETURN 
1760 DEFFN'04! SELECT PRINT 005! PRINT HEX (03) ! PLOT <1,,C>,<,,S> 
1770 PRINT “POSITION LABEL WHERE IT WILL NOT INTERFERE WITH GRAP 
H" 

1780 P$~" " ! KEY IN P$, 1800, 1810 
1790 GOTO 1780 

1800 IF P*-HEX(OD> THEN 1880! IF P$-HEX(02) THEN 1830! IF P$-HEX<0 
8) THEN 1840: PLOT <, , P$>, <13, , LJ>! GOTO 1780 

1810 IF P$-HEX ( 05 ) THEN 1850; IF P$-HEX(06) THEN I860; GOTO 1870 

1 820 PLOT <0 , -20 , U> , <-999 , 0 , U> ! GOTO 1 780 

1830 PLOT <13, 0,U>! GOTO 1780 

1840 PLOT <- 1 3,0, U> ! GOTO 1780 

1 850 PLOT <0 , 20 , U> ! GOTO 1780 

1860 PLOT <0,-20,U>! GOTO 1780 

1870 STOP 

1880 RETURN 

1 8 9 0 P L 0 T < 1 , , C > , < 1 4 , , 8 > , < , , E $ > 

1900 PLOT <-238, -25, U>! PLOT <,,J$> 

1 9 1 0 P LOT < - 140,-2 5 , U > ! PL 0 T < , , D $ > 

1920 PLOT <-98, -25, U>! GQSIJB 1730! PLOT <18,,G$> 

1930 PLOT <-214,-25, U>! GOSUB 1740! PLOT <18, s , HO 
1940 PLOT <-228, -25, U>! GOSIJB 1750! PLOT <18, , I$>, <, , R> 

1950 RETURN 

1960 PRINT “LOAD DATA TAPE IN UNIT 10A" ! INPUT “FILE NUMBER TO BE 
LOADED", Ni 

1970 SELECT TAPE 10A 

1980 P9-0! REWIND ! FOR K*1 TO 50! P9-P9+1 ! NEXT K! IF Nl=l THEN 1990 
! N-Nl-l ! SKIP NF 

1990 DATA LOAD " LTEST" ! DATA LOAD D$ , J$ , Z 1 , D ( 8 ) , Q ( ) , H ( ) , P ( ) , Q 1 ( ) ! 
IF END THEN 2000; GOTO 2000 
2000 SELECT PRINT 215(132) 

2010 PRINT “RSECS PPC02 DATA: FILE NUMBER" ; Nl; PRINT 
2020 PRINT “TEST DATE-" ; DY> 
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2030 PRINT "TYPE OF CASE» M JJ* 

2040 PRINT “VALUE OF X AT THE ORIGIN-" ) 21 

20 SO PRINT "TOTAL NUMBER OF DATA POINTS-" ; D<8) ! PRINT 

2060 PRINT y " TIME" r " ARSIN" y " UPCAN" , " LOWCAN" 

2070 PRINT " I V' Q < I > " y " H ( I ) " y " P ( I > " , " Q1 < I > " 

2080 FOR 1=1 TO DCS) 

2090 PRINT I,QCI),H(I)yPCI)yQlCI) 

2100 NEXT I 

2110 SELECT PRINT 00SC64) 

2120 STOP 


LIOHCRT COMPUTER PROGRAM 


File Name "LIOHCRT" 

Abstract "LIOHCRT" plots the instantaneous removal rate of CO 2 

from the LiOH canisters in lbm/min. The program pro- 
duces a separate graph for both the upper and lower 
LiOH canisters using the LiOH test data previously 
saved on cassette tapes by the "LIOHPPC" program. 

Figures 7 and 8 provide examples of the plots produced 
by the program. The program is designed for use with 
the Wang 2200 series computer system with a flatbed 
plotter. 

Program Description 

The program takes data previously saved on cassette tapes by program 
"LIOHPPC" and calculates a CO 2 removal rate from the LiOH canisters. 

The formula used in the program is as follows: 

M C 0 2 (REMOVAL) = 0.07931* V *A P/(T+459.67) 

Withv in CFM, p in mmHg, and T in deg F, the formula calculates the 
instantaneous removal rate of CO 2 in lbm/min. A full discussion of 
the derivation of this formula is included in Appendix 1. 

The upper and lower LiOH canisters were changed out alternately with 
unique temperature and flow rates associated with each canister. The 
temperatures and flow rates used were the averaged values taken at 
canister insertion and removal, and are available in Memorandum EC2- 
.77-185, "RSECS LiOH Test Data", dated 11-8-77. The length of time 
for a LiOH canister may encompass more than one file on the data tapes. 
The program first plots the upper canister removal rate and then the 
lower canister removal rate. Therefore, after loading a data tape 
in the tape drive unit, the program requests a file number, first and 
last data points to be plotted, and the temperature and flow rate 
associated with the first upper LiOH canister. The plotting is per- 
formed incrementally as data points are specified. 

Any one file may not include all the data points associated with a 
canister. Therefore, another file, the first and last data points 
to be plotted and the temperature and flow rate must again be specified 
until all data points associated with the canister are reduced and 
plotted. When a specific canister is completed, either the next set 
of data points in a file or another file is required for the next can- 
ister. Thus, the file number, first and last data points to be plotted, 
and the temperature and flow rate for the new canister is requested by 
the program. In this way, the user will plot CO 2 removal rates for 
all of the upper canisters in a specific case. 

The program provides labels for the plots, but requires that the user 
position the label where it will not interfere with the graph. Next, 
the program requests that the user replace the graph paper on the 
flatbed plotter before plotting the lower LiOH canister. After the 
graph paper has been changed, the procedure outlined for the upper 


canister is repeated for the lower canister. Figures 7 and 8 exemplify 
plots produced by the program for the 7 man case. 

The program interfaces the input data with the generalized plot routine 
as a subroutine in order to produce the plots. A listing of the program 
"LIOHCRT" is. provided in Table VII. 
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TABLE VII LIOHCRT PROGRAM LISTING 


10 REM L I OHCRT- PROG RAM PLOTS COE REMOVAL RATES FROM LI OH BEDS 
BASED UPON DATA FROM TAPES GENERATED BY LIOHPPC PROGRAM 
20 COM X9 ( 1 50 ) , Y9 ( 1 SO ) „ C ( 1 0 ) , X$25 , Y$S5 , P$ 1 , C$30 , D$EG 
30 COM E$20,F$26,,Z1,G$25,H$25, I*, J$,R*20 
40 COM Q( 150), H( 150), P ( 150), ,D( 10), 01(150) 

50 E$=" LI OH TEST RESULTS" 

00 PRINT "LOAD DATA TAPE ON UNIT B "’.INPUT "FILE NUMBER TO BE. LO 
ADED " „ N 1 "SELECT TAPE 10B : P9=G: REWIND l FOR K*1 TO 50:P9=P9+1: 
NEXT K:IF Nl=l THEN 70: N=N1-1 : SKIP NF 
70 DATA LOAD "LTEST" : DATA LOAD D$, J$, Z 1 , D(S ) , Q ( ) , H( ) , P ( ) , Q1 ( ) : 

IF END THEN SO: GOTO 30 

SO SELECT TAPE 10A: K6=D ( 8) : X$= " TIME IN MINUTES" : Y$=" COE REMOVAL 
RATE IN LB /MIN" 

90 r$~d$; X 0=400: Y0= . 005: x l * l : yi= l : c 1 *z 1 : C2=0: si =z l : S2=z l +3600: Tl 

=0: T2=. 045: LI *2: L2*l : GOSUB ' 00 

100 PRINT "PLOTTING UPPPER CANNISTER " 

110 G$=" UP CAN COE REMOVAL RATE" 

120 INPUT "FIRST AND LAST DATA POINTS*" , II , 12 
130 INPUT "TEMP AND FLOW*" , T, V5:K6=I2 
140 Q$= "FIGURE 1" 

150 FOR 1=11 TO IE: X9( I ) =Q ( I ) : Y9 ( I ) =V5#. 07931 / (459. 67+T )*(H(I)~P 
(I)):NEXT I 

1 60 P*=HEX ( 00 ) : D= 1 : K= 1 : K 5=K6 : LS= 1 1 : L9= I s : GOSUB ' 03 
170 FOR 1=1 TO 1 50: X9( I ) =0: Y9( I ) =0: NEXT I 

ISO INPUT "MORE DATA POINTS? (0=N0, 1=YES) " , P3: IF P3=l THEN ISO 
190 INPUT "ANOTHER FILE? (0=NQ, 1=YES) " , A; IF A=0 THEN 200: GOSUB 
1790*. GOTO 120 

EOO INPUT "DO YOU WISH TO LABEL GRAPH? (0=NO, 1=YES)",P:IF P=Q T 
HEN 210: GOSUB "04: GOSUB 17S0 

S10 PRINT "MOW PLOTTING LOWER CAN GRAPH, CHANGE PAPER": STOP : GOS 
UB '00 

250 INPUT "FIRST AND LAST DATA POINTS*" , II, IE 

230 INPUT "TEMP AND FLOW*" , T , V5; K6=IE: Gk ~ "LOW CAN COE REMOVAL RA 
TE" 

240 FOR 1=11 TO IE:X9<I)=Q(I):Y9(I>=V5*. 07931/(459. 67+T)*(H< I )-Q 

1(D): NEXT I 

250 Q$=" FIGURE 2“ 

EGO P $ =HE X ( 0 1 ) : D = 1 : K = 1 : K 5 =K 6 : LS = 1 1 : L9 = 1 2 : GOSUB '03 
E70 FOR 1=1 TO 150:X9(I)=0:Y9(I)=0:NEXT I 

ESO INPUT "MORE DATA POINTS? ( 0=N0 , 1 = YES ) " , P 3: IF P3=l THEN EEO 
290 INPUT "ANOTHER FILE? (G=NQ, 1=YES)",A:IF A=0 THEN 300: GOSUB 
1790: GOTO EEO 

300 INPUT "DO YOU WISH TO LABEL GRAPH? (0=N0, 1-YES) ", P: IF P=0 T 
HEN 3S0 

310 GOSUB '04: GOSUB 17E0 
3E0 STOP 

330 DEFFN ' 00”. Q=0: GOSUB 340: PLOT <, ,R>: RETURN 
340 SELECT PRINT 005: PRINT HEX (03 )*. PRINT SPRINT *. 
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PRINT "#*# GENERATING PLOTS OF " ;E$;“ ***" 

350 PRINT ; PRINT : PRINT " IS PLOTTER POWER ON ?" SPRINT " IS PLOT 

PAPER LOADED ? " i PRINT " IS CHART SWITCH ON HOLD ? PRINT 11 IS 

PEN IN DOWN POSITION AND CAP OFF ?" 

360 PRINT " ARE SCALE LIMITS SET <0 y 0 BY 10,10) ?" 

370 INPUT " IS ANSWER TO ALL ABODE QUESTIONS YES ? ( YES»I, N0=0> " 


Q! 


IF QOl THEN 340 


• 5> 9 w 


380 SELECT PLOT 414! PLOT <1, r C>,<12, , S>! Q=0 
390 F1-100/X0 ! F2-100/Y0 : GOSUB 630 ! PLOT <1, y C>, 

400 RETURN 

410 DEFFN / 03: SELECT PRINT 005! PRINT HEX C 03 > ! 

PLOT < y y R>, <100«Xl y 100*Y1 y U> y <1 , , C> 

420 X4 f Y4, E ? E3 y E4 y E6 y E8 y E7--0! Ul-l; U3-0 
430 K~1 ! X4 y Y4-0 
440 Kl-KlW-1 

450 U2--0 ! X=X9 < K ) ! I F X<Z1 THEN 500; IF X>S2 THEN 500; Y=Y9(K) ! IF Y> 
T2 THEN 5001 GOSUB '125; IF U2=l THEN 520; X-X9 ( K ) -Cl ; Y-Y? (1-0 --C2: X4 
. Y4-0 

460 IF K>K1 THEN 480: PLOT < y f R>. <iOO*Xl y 100*Y1 , U> y <F1*X, F2*Y y U> 

4 7 0 P L 0 T , D > : G 0 T 0 490 

480 GOSUB ' 122 <Xj Y y X4 y Y4 ) 

490 PLOT <-X*Fi y ~Y*F2 y U>; GOTO 520 
500 IF K/6-W THEN 510! GOTO 520 
510 W*W+1 

520 K«K+l: IF IO--K5+1 THEN 530; GOTO 450 

530 X4 y Y4, E, E3, E4 y E6 y E8* E7-0: Ul»2! PLOT < y y R> y <100#X1, 100*Y1, U> 

540 FOR I-L8 TO L9 

550 IF X9 ( I ) <C1 THEN 600 

560 X=X9 ( I ) -Cl ! Y~Y9 ( 1 ) -C2 

570 IF X>3600 THEN 600: IF Y>T2 THEN 600 

580 IF lOLS THEN 590 i U 1 ~ 1 : GOSUB ' 122 (X, Y,X4, Y4 ) 1 U1 «2 : GOTO 610 
590 GOSUB 1 22 ( X y Y y X4 , Y4 ) : GOT 0 610 
600 U3«l 
6.10 NEXT I 

620 PLOT < y , U>: PLOT < y y R>: RETURN 

630 Al-Fl W ABS ( Sl-Cl ) ; A2«F1*ABS<S2~C1 ) ; B1«F2#ABS < T1-C2 ) ; B2-F2#ABS 
< T2-C2 ) 

640 PLOT < y , R> y <100*' XI , 100*Y1 y U>, <-Al y 0 y U>, <A1+A2 y 0, D>, <-A2 y ~Bl y 
U> y <O y Bl+B2 y D> 

650 M 5 - ( A B S ( S 1 •••• C 1 ) +ABS CS2-C1 ) >/X0: NS-«ABS(T1-C2)+ABS(T2-C2) )/Y0 

660 K-0 

670 S3-S1-X0 

680. PLOT <-~3 y •••• ( B.1. + B2 ) y U> 

690 N5-INT <NE*. 5 ) 

700 FOR 13-1 TO N.5+1 
710 PLOT <6 r 0,D>,<, y U> 

720 IF I3-N5+1 THEN 730: PLOT <~6 y 100 y U> 

730 NEXT 13 
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740 PLOT <-- ( A 1 + 3 ) , - ( B2+6 ) , U> 

750 MS«INT(M5+. 5) 

760 FOR I 4~ l TO M5+1 
770 PLOT <0,12,0>,<,,U> 

780 IF I4”M5+ 1 THEN 790; PLOT <100,-12,U> 

790 NEXT 14 

800 IF L_l=2 THEN 810! PLOT <--( A1+A2+24 ) , 20 y U>! GOTO 820 
810 PLOT <-<Al+A2+24),-32,U> 

820 FOR 1=1 TO M5+1 
830 IF I>M5+1 THEN 940 
840 S3=S3+X0 

850 IF ABS ( S3 ) >-= 1000. THEN 870; IF ABS ( S3 ) >=100 « THEN 880; IF ABS< S 

3 ) >=10. THEN 890; IF ABS ( S3 >> = 1 » THEN 900 

860 CONVERT S3 TO S3*, <-.###> I GOTO 910 

870 CONVERT S3 TO S3*, <-####>: GOTO 910 

880 CONVERT S3 TO S3$, ( ) I GOTO 910 

890 CONVERT S3 TO S3* ,<-##.#>! GOTO 910 

900 CONVERT S3 TO S3*, <-#.##)! GOTO 910 

910 IF KOO THEN 920: PLOT <, , S3*> ! GOTO 940 

920 IF I-IOINT < A1 + . 5 ) /100 THEN 930; PLOT <100, 0, U> I GOTO 940 

930 PLOT <100-60, 0, U>, <, , S3*> 

940 K=K+1 : NEXT I 

950 IF L2«l THEN 960; PLOT <~< A2+18 ) , 0, U> : GOTO 970 
960 PLOT <-<A2+95>,0,U> 

970 IF Ll- : 1 THEN 980: PLOT < 0, - < Bl-26 ) , U> ! GOTO 990 
980 PLOT <0, ~<Bl+26) y U> 

990 K®0 

1000 T3=T1--Y0 

1010 FOR 12=1 TO NS+1 

1020 IF I2>N5+1 THEN 1130 

1030 T3--T3 + Y 0 

1040 IF ABS ( T3 ) >- : 1 000 « THEN 1060: IF ABS < T3 ) >=.100» THEN 1070; IF ABS 

( T3 ) >= 1 0 - THEN 1 080 : I F ABS ( T3 ) >« 1 » THEN 1 090 

1050 CONVERT T3 TO T3*, <-.##»)! GOTO 1100 

1060 CONVERT T3 TO T3*, <-####): GOTO 1100 

1070 CONVERT T3 TO T3*, (-##.#.■>: GOTO 1100 

1080 CONVERT T3 TO T3*, <-H#. # > : GOTO 1100 

1090 CONVERT T3 TO T3*, I GOTO 1100 

1100 IF KOO THEN 1110: PLOT <, , T3*> ! GOTO 1130 

1110 IF 1 2 -- 1 <> I N T < B 1 + « 5 ) / 1 0 0 THEN 1120: PLOT <0,100,U>! GOTO 1130 
1120 PLOT <“60, 100, U>, <, ,T3*> 

1130 K«K+1!NEXT 12 

1.140 IF X$=" " THEN 1200: PLOT <, , R>, <100*X1-A1, 100*Y1, U> 

1150 IF Ll*2 THEN 1160; PLOT <0, SO, U> I GOTO 1170 
1160 PLOT <0,-60-, U> 

1170 IF A 200 THEN 1180! PLOT <Al/5, 0, U> I GOTO 1190 
1180 PLOT <Al+A2/2« 25,0, U> 

1190 PLOT <, , X*> 
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1200 IF Y$-" 11 THEN 1270 

1210 PLOT <, , R>, <100*X1 , 100*Y1-B1 , U> 

1220 IF L.2-1 THEN 1230: PLOT <90, '0, U>! GOTO 1240 
1230 PLOT <-90, 0, U> 

1240 IF B2<>0 THEN 1250; PLOT <0, 2*01/3, U>; GOTO 1260 
1 2 SO PLOT <0 , B 1 + B2# 10.5/12, U> 

1260 P LOT <0,-20, S > , < , , Y $ > , < 13,, S > , < . , R > 

1270 RETURN 

1280 DEFFN ' 122 < U, 0, X4, Y4 ) 

1290 X3-U!Y3*V 

1300 D1»X3~X4: D2= ; Y3-Y4: X4=X4+D1 J Y4=Y4+D2 
1310 IF Ul-2 THEN 1370 

1320 IF U2-1 THEN 1350: PLOT <F1*D1 , F2*D2, U>, <, , D> 
1330 PLOT <,„[)>: PLOT <, , U> ! GOTO 1360 
1340 GOSUB 1550: W-W+l : PLOT <,, U>;GOTO 1360 
1350 PLOT <F1*D1 , U3, U>, <, , U>, <, , >, <, , U>: 

E3-.TNT (U3 ) : E7-INT ( F1*D 1 ) : GOTO 1500 
1360 E3-INT ( F2*D2 ) ! E7-INT ( F1*D1 > I RETURN 
1370 E1=F2#Y.3~INT (F2*Y3 ) 

1380 E-E1+E4 


1390 P8=F2#D2+E: P9=INT(P8) 

1400 E5=F1*X3-INT ( F1*X3 ) 

1410 E6=ES+E8 

1420 S8-F1*D1+E.6 ! S9~INT(S8) 

1430 IF U2-0 THEN 1440: PLOT <, , U>, <89, 0, U>: GOTO 1460 

1440 IF U3=l THEN 1450; PLOT <, , D>„ <89 , P9, D> : GOTO 1460 

1 450 PLOT <, , LJ> , <89 , P9 , U> f U3=0 

1460 IF IJ2-1 THEN 1470; E3-E3+P9 

1470 E4-F2*Y3~E3 

1480 E7-S9+E 7 

1490 E8~F1*X3'*E7 

1500 RETURN 

1510 DEFFN / 125 ; IF ( Y-C2 ) /Y0 + YK9 . 999THEN 1S20:U2=1 
1520 IF ( Y-C2 ) / Y0+Y1>. 001 THEN 1530: U2-1 
1530 IF (X--C1 )/X0+Xl<9„ 999THEN 1S40:U2=1 
1540 RETURN 

1550 ON <VAL<P*)+1>G0SUB 1560, 1570, 1580: RETURN 

1560 PLOT <0, 8, U>, <-7, -12, D>, <14, 0, D>, <~7, 12, D>, <0, -8, U > \ 


RETURN 


1570 PLOT <0, 6, U>, <-4, -12, D>, <10, 8, D>, <-12, 0, D>, <10, -8, D>, <-4, 12 
, D>, <0,-6, U>! RETURN 

1580 PLOT <-4, 7, D>, <8, 0, D>, <-8, -14, D>, <8, 0, D>, <-4, 7, D>, <-7, 4, D>, 
<0,-8, D>, <14,8, D>, <0, -8, D>, <-7,4, D>, <0, 0, U>! RETURN 
1590 DEFFN / 04! SELECT PRINT 005: PRINT HEX ( 03 >} PLOT <1,,C>,<, ,S> 
1600 PRINT "POSITION LABEL WHERE IT WILL NOT INTERFERE WITH GRAP 
H" 


1610 P$~" " ; KEY IN P$, 1630, 1640 
1620 GOTO 1610 
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1630 IF r$~HEX(0D) THEN 1710! IF P$~HEX(02) THEN 1660! IF P$=HEX(0 
8 ) THEN 1 670 ! PLOT < , , P$> , < 13 , , U> ! GOTO 1610 

1640 IF P$~HEX ( 05 ) THEN 1680! IF P$»HEX< 06) THEN 1690! GOTO 1700 

16 SO PLOT < 0 y - 2 0 , U > , < -- 9 9 9 , 0 , U > ! GOT 0 1610 

1660 PLOT <13, 0, U>! GOTO 1610 

1 670 PLOT <-- 13,0, U> ! GOTO 1610 

1680 PLOT <0,20,U>: GOTO 1610 

1690 PLOT <0,-20,U>! GOTO 1610 

1700 STOP 

1710 RETURN 


1720 
1730 
1740 
17 SO 
1760 
1770 
1780 


PLOT 

PLOT 

PLOT 

PLOT- 

PLOT 

PLOT 


f ' 

, , W 

,•=•.14, , 

O , 

Q$> 

112, 

~2S, U> 

: PLOT 

<,,E* 

238, 

~2S, LO 

! PLOT 

<f f TT 

140, 

-2S, U> 

: PLOT 

<, , R$ 

98, ••• 

2S, u>: 

PLOT ■: 

, g$> 


RETURN 


1790 PRINT “LOAD DATA TAPE ON UNIT B 11 ! INPUT "FILE NUMBER TO BE 
LOADED", N1 ! SELECT TAPE 10B! P9--0! REWIND ! FOR K=1 TO 50;P9=P9 + l! 

NEXT K‘. IF Nl»l THEN 70! N*Nl-l! SKIP NF 
1800 DATA LOAD "LTEST" ! DATA LOAD D$, J», Z1 , D < 8 ) , Q ( ) , H O , P < > , 01 ( > : 


IF END THEN 1810! GOTO 1810 


1810 SELECT TAPE 10A ! K6«D ( 8 ) ! X*«" TIME IN MINUTES" ! Y$="C02 REMOOA 
L RATE IN LB /HR" 

1820 RETURN 


ow 

rtf 


■sssss* : 
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LIOHMAC COMPUTER PROGRAM 


File Name "LIOHMAC" 

Abstract "LIOHMAC" plots the accumulation of mass of C02 in the 
— lioh canisters in lbm. The program produces a separate 

plot for the upper and the lower canisters using the Li OH 
test data previously stored on cassette tapes by the 
"LIOHPPC" program. Figures 9 and 10 provide examples of 
the plots produced by the program. The program is de- 
signed for use with the Wang 2200 series computer system 
with a flatbed plotter. 

Program Description 

The program utilizes data previously saved on cassette tapes by program 
"LIOHPPC" to calculate an accumulation of mass of CO2 in the LiOH can- 
isters. The basic formula used in the program is as follows: 

MASS ACCUMULATED (lbm) = RATE OF MASS REMOVAL (lbm/mln)* 

DELTA TIME (min) 

RATE OF MASS REMOVAL (lbm/min)= 0 . 07931*V*A ?/ (T+459 . 67) 

where V is in CFM, P, is in mmHg, and T is in degrees F. The formula for 
rate of removal is the same as used in program "LIOHCRT" and is fully 
discussed in Appendix 1. The program "LIOHMAC" functions exactly the 
same way as "LIOHCRT" with the expanded formula being the primary dif- 
ference. Therefore, the user may refer to the program description of 
"LIOHCRT" for an explanation of the required input data and program 
operation. 

Figures 9 and 10 exemplify plots produced by the program for the 7 man 
case. The program interfaces the input data with the generalized plot 
routine as a subroutine in order to produce the plots. A listing of 
the program "LIOHMAC" is provided in Table VIII. 


original ways 
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TABLE VIII LIOHMAC PROGRAM LISTING 




10 REM LIOHMAC -PROGRAM PLOTS 002 ACCUMULATION IN LIOH BEDS 
BASED UPON DATA FROM TAPES GENERATED BY LIOHPPC PROGRAM 
20 COM X9 ( ISO > . Y9 < ISO ) , C < 10 ) , X$25. Y$25,P$1 , C$30, D$20 
30 COM E$20, F$20, Z1 , G$25, H$2S, 1$, J$„ R$20 
40 COM Q < 1 SO > , H (ISO), P ( 1 SO ) , D( 10) , Q1 ( ISO) 

SO E$-' : " LIOH TEST RESULTS” 

60 PRINT "LOAD DATA TAPE ON UNIT B ” 1 INPUT "FILE NUMBER TO BE L.0 
ADED" , N1 ; SELECT TAPE 10B; P9 : ~0; REWIND : FOR K~1 TO 50! P9=P9U ; 

NEXT K : IF N 1 ™ 1 THEN 70 ! N™N1™1 ! SKIP NF 
70 DATA LOAD " L..TEST" ! DATA LOAD D$. J$, Z1 , D < 8 ) , Q O , H < ) , P ( > , Q1 C ) ! 

IF END THEN SO! GOTO 80 

80 SELECT TAPE 1 0 A ! K6™D ( 8 ) ! X$~ " TIME IN MINUTES" ! Y$™"C02 ACCUMULA 
TION IN LBS" 

90 R$™D$ ! X0-400! YO™. 5; Xl = l! Yl = l! C1”Z1 ! C2^0 ! S1“Z1 ! S2“Zl + 3600 ! T1=0 

! T2™4 u S ! L_ 1 =-2 ! L2-1 ; GOSUB ' 00 

100 PRINT "PLOTTING UPPPER CANNISTER" 

110 G$™ " UP CAN C02 ACCUMULATED MASS" 

120 INPUT "FIRST AND LAST DATA POINTS™" , II, 12 
130 INPUT "TEMP AND FLOW--" , T, MS; K6™12 
140 Q$-=" FIGURE 1" 

ISO FOR I ~ 1 1 TO 12! Y9(Ii™l)="Y9(Il™l)+M 

1 60 X9 ( I ) -Q ( I ) ! Y9 ( I ) :==V5* . 0793 1 / < 4S9 . 67+T > * ( l-f < I ) -P ( I ) ) * ( ( ( Q ( I ) +Q ( 
I 1 ) ) /2 ) - ( ( Q ( I •- 1 > + Q ( I ) ) /2 ) ) +Y9 (1-1 ) ! NEXT I 
170 M--Y9 ( 12 ) 

1 80 P$= : HEX (00 ) ! D~ 1 ! K- 1 ! KS™-K6 ! L8™ 1 1 ! L9™ 1 2 ! GOSUB y 03 
190 FOR I~1 TO 1 SO ! X9 ( 1 ) ™() ! Y9 ( I ) ™0 ! NEXT I 

200 INPUT "ANOTHER FILE? (0*N0, 1™YES ) " , A ! IF A=0 THEN 210; GOSUB 
I860! GOTO 120 

210 M™0: INPUT “MORE DATA POINTS? < O-NO, 1™YES ) " , P3 ! IF P3=l THEN 1 
20 

220 INPUT “DO YOU WISH TO LABEL GRAPH? (0™NO, 1~YES) " , P! IF P™0 T 
HEN 230! GOSUB '04: GOSUB 1790 

230 PRINT "NOW PLOTTING LOWER CAN GRAPH, CHANGE PAPER" ! STOP I GOS 
UB '00 

240 PRINT "LOAD DATA TAPE ON UNIT B "I INPUT "FILE NUMBER TO BE L 
OADED" . N1 : SELECT TAPE 10A : P9=0 1 REWIND IFOR H>1 TO SO : P9-P9+1 : 
NEXT K: IF N 1 =- 1 THEN 70 1 N™N1-1 ; SKIP NF 
250 DATA LOAD " L.TEST" : DATA LOAD D$ , J$ , Z1 , D ( 8 ) , Q ( ) , H ( ) , P ( ) , Q 1 ( ) ; 

IF END THEN 260: GOTO 260 

260 SELECT TAPE 10A I K6~D ( 8 ) : X$~" TIME IN MINUTES" i Y$“" C02 ACCUMUL 
AT I ON IN LBS" 

270 INPUT "FIRST AND LAST DATA POINTS™" , II » 12 

280 INPUT "TEMP AND FLOW™" , T, US: K6™12: G$“" LOW CAN C02 ACCUMULATE 
D MASS" 

290 FOR I»Ii TO 1 2 : Y9 ( 1 1 - 1 ) = : Y9 ( 1 1 - 1 ) +M 

300 X9( I ) --0 ( I ) ; Y9( I )™V5». 07931/(459. 67+T)*(H( I ) ~ P ( I ) ) * ( ( <Q( I )+Q( 
I 1 ) ) /2 ) ( ( « ( I ™ 1 )+Q( I ) ) /2) )+Y9( I- l ) ; NEXT I 
310 M--Y9 (12) 

320 Q$--" FIGURE 2" 
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3 3 0 P $ a H E X ( 0 i } ! D " 1 : K » 1 ! K S K 6 : !... 8 1 1 : L 9 * 1 2 ; 6 0 S U B / 0 3 
34 0 FOR I 1 TO 1 SO ! X9 < I > «0 : Y9 ( I ) «0 1 NEXT I 

3 5 0 I N P U T " A N 0 T H E R F I L. E ? < 0 « N 0 , 1 “ Y E S ) % A : I F A ~ 0 THEN 3 6 0 ! G 0 S U & 
•I860: GOTO 270 

360 M=~0: INPUT "MORE DATA POINTS? <0*N0, 1~YES> " , P3i IF P3=l THEN 2 
70 

370 INPUT "DO YOU WISH TO LABEL GRAPH? <0«ND, I -YES) " y P! IF P~0 T 
HEN 390 

380 GOSUB 7 04 i GO SUB 1790 
390 STOP 

400 DEFFN 7 00: 0-0 : GOSUB 410: PLOT <,,R>: RETURN 
410 SELECT PRINT 005! PRINT HEX <03 ): PRINT ! PRINT : 

PRINT “*** GENERATING PLOTS OF " ; E$J " #*#" 

420 PRINT : PRINT ! PRINT " IS PLOTTER POWER ON ?" SPRINT " IS PLOT 
PAPER LOADED ? "{PRINT " IS CHART SWITCH ON HOLD ? "SPRINT " IS 
PEN l'N DOWN POSITION AND CAP OFF ?" 

430 PRINT " ARE SCALE LIMITS SET (0,0 BY 10,10) ?" 

440 INPUT " IS ANSWER TO ALL ABOVE QUESTIONS YES ? <YES»1,ND*0>" 
,0; IF 0O1 THEN 410 

4 SO SELECT PLOT 414: PLOT <1 , , O, <12, , S> ! 0=0 

460 Fl*100/X0: FS-IOO/YO: GOSUB 700: PLOT <1, , C>, <, , SI- 

470 RETURN 

480 DEFFN 7 03! SELECT PRINT 005: PRINT HEX (03): 

PLOT < g „ R>, <100*X1, 100*Y1,U>, <1, , C> 

490 X4„ Y4, E, E3, E4 , E6, EB, E7=0 1 U1---1 : U3=0 
SCO K™ 1 1 X4 , Y4~0 
S10 K.1. -•• !•<: W -“l 

520 U2«0:X«X9<K): IF X<Z1 THEN 570: IF X>S2 THEN 570 : Y~Y9 ( K ) : IF Y> 
T2 THEN 570; GOSUB ''125: IF U2=l THEN 590 ; X=X9 ( K ) -Cl : Y=Y9 ( K ) -C2 : X4 
, Y4~0 

530 IF IOK1 THEN 550: PLOT <, , R>, <100*X1, 100*Y1 , U>, <F1*X, F2#Y, U> 

5 4 0 P L 0 T < . y D > G 0 T 0 560 
550 GOSUB ' 122 ( X, Y, X4 y Y4 ) 

560 PLOT <~X*F1,~Y*F2,U>! GOTO 590 
570 IF K/5=W THEN 580: GOTO 590 
580 W®W+1 

590 K-K+l : IF K>--"K5+1 THEN 600; GOTO 520 

600 X4, Y4 , Ej E3, E4,E6,E8, E7~0 : U1--2 : PLOT <, ,«>, <1(>0*X1 , 100#Y1, U> 

610 FOR I«L8 TO L9 

620 IF X9 ( I )<C1 THEN 670 

630 X~X9 ( I > -C 1 ! Y“Y9 ( I ) ~C2 

640 IF X>3600 THEN 670; IF Y>T2 THEN 670 

650 I F I OL8 THEN 660 ; U 1 « 1 ! GOSUB 7 1 22 (X,Y, X4 , Y4 ) ! U1 ™2 : GOTO 680 
660 GOSUB 7 122 (X, Y, X4 y Y4 ) : GOTO 680 
670 U3"“l 
680 NEXT I 

6 9 0 P L 0 T < , , U > : P L 0 T < , , R> I R E T U R N 

700 A1=F1*ABS(S1~C1 ) : A2«Fi*ABS(S2~Cl ) ! B1-F2*ABS ( T1-C2 ) l B2=F2*ABS 
( T2-C2 ) 


TABLE VIII CONT 


710 PLOT <. , R > , <1 00#X1, 100*Y1 , U>, <~A1 , 0, U>, CA1+A2, 0, D>, <~A2, -B1 , 
U>.<0,Bi+B2,D> 

720 M5= < ABB ( SI -Cl > +ABS ( S2-C1 ) ) /XO ! NS- ( ABB ( T1-C2 ) +ABS < T2-C2 ) > /YO 


730 K-0 
740 S3~81~X0 

7 SO PLOT <~‘3 , - < B 1 +B2 ) , U> 

760 NS- I NT CNS+. S) 

770 FOR 13-1 TO NS+1 
7 8 0 P L. 0 T < 6 . 0, D > , < y , U > 

790 IF I 3- NS+1 THEN 800: PLOT <-6,100,U> 
800 NEXT 13 


8 1 0 PLOT <- ( A 1 +3 ) , - ( B2+6 > , U> 
820 MS- IN) ( MS+ . S) 

830 FOR 14-1 TO MS+1 


840 PLOT <0, 12, D>, <, , U> 

8S0 IF I 4 -MS+1 THEN 860: PLOT <100, -12, U> 

860 NEXT 14 

870 IF L.l-2 THEN 880: PLOT <- ( A1+A2 + 24 ) , 20, U> ! GOTO 890 
880 PLOT < - ( A 1 + A 2 +• 2 4 > ,-32 , U > 


890 FOR 1-1 TO MS+1 
900 IF I>MS+1 THEN 1010 
910 S3-S3+X0 

920 IF ABB < S3 ) >- 1 000 . THEN 940: IF ABB ( S3 > >- 1 00 THEN 950: IF ABS< S 

3 ) >== 1 0 . THEN 960 : I F ABB C S3 ) >- 1 » THEN 970 

930 CONVERT S3 TO S3$, (-.««#) : GOTO 980 

940 CONVERT S3 TO S3*. (-•###«>; GOTO 980 

9S0 CONVERT S3 TO S3*. (-###.): GOTO 980 

960 CONVERT S3 TO S3*. <-##.#)! GOTO 980 

970 CONVERT S3 TO S3*. <-«.##): GOTO 980 

980 IF KOO THEN 990: PLOT <, f S3*>:G0T0 1010 

990 IF 3! ~ 1 < > I NT < A 1 + » S ) / 1 0 0 THEN 1000: PLOT <100, 0, U>: GOTO 1010 


'* f y v O 


3*> 


1000 PLOT <100-60, 0,U> 

1010 K-K+-1 : NEXT I 
1020 IF L2=l THEN 1030! PLOT 
1 030 PLOT <- < A2+9S ) , 0 , U> 

1040 IF Ll-1 THEN 10 SO: PLOT < 0 
1 0 SO PL 0 T <0 , - < B 1 + 26 ) , U> 

.1060 K-0 

1070 T3-T1- YO 

1080 FOR 12-1 TO NS+1 

1090 IF I2>NS+1 THEN 1200 

1100 T3-T3+Y0 

1110 IF A B S ( T 3 ) > - 1 0 0 0 » T HE N 1130 
( T3 ) >-10 « THEN 11 SO: IF ABS(T3>>- 
1 1 2 0 C 0 N V E R T T 3 T 0 T 3 *,(-.# # # ) : 

1130 CONVERT T3 TO T3*, <-###«): 

1140 CONVERT T3 TO T3*, <-##«.) ! 

1150 CONVERT T3 TO T3*, 

1160 CONVERT T3 TO T3*, (-#.«#): 


< A2 + 18 ) , 0, U> 


< B 1 -26 > , U 


I IF ABSCT3) 
1 » THEN 1160 
GOTO 1170 
GOTO 1.170 
GOTO 1170 
GOTO 1170 
GOTO 1170 


GOTO 1040 
IGOTO 1060 


-100. THEN 1140: IF ABS 
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1170 IF l«>0 THEN 1180! PLOT <, , T3*>5 GOTO 1200 

1180 IF 1 2 - 1 <> I N T ( B 1 >, 5 ) / 1 0 0 THEN 1190: PLOT <0, 100, U>: GOTO 1200 
1 1 90 PLOT <-60, 1 00 , U> , < , y T3$> 

1200 K=K+1 ! NEXT 12 

1210 IF Xf*»" " THEN 1270! PLOT <, , R>, <100*X1-A1 , 100*Y.t, U> 

1220 IF L_l”2 THEN 1230; PLOT <0, 50. U>! GOTO 1240 
1230 PLOT <0, -60, U> 

12.40 IF A2O0 THEN 1250; PLOT <A1/S, 0, U>; GOTO 1260 
1250 PLOT <A1+A2/2.2S,0,U> 

1260 PLOT <, ,X$> 

1270 IF Y$~" 11 THEN 1340 

1280 PLOT <, . R> , <100*X1 , 100«Y1-B1 , U> 

1290 IF L2»l THEN 1300: PLOT <90, 0. U>: GOTO 1310 
1300 PLOT <-90,0. U> 

1310 IF B2O0 THEN 13201 PLOT <0,2#Bl/3,(J>i GOTO 1330 
1320 PLOT <0,81+62*10. 5/12, U> 

1330 PLOT <0, -20, S>, <, , Y$>, <13, , S>, <, , R> 

1340 RETURN 

1350 DEFFN ' 1 22 ( U , M , X4 , Y4 ) 

1360 X3~(J; Y3™V 

1370 D1-X3-X4 ; [)2™Y3“Y4: X4-X4+D1 1 Y4-Y4+D2 
1380 IF U 1 “2 THEN 1440 

1390 IF U2~l THEN 1420: PLOT <F1*D1 , F2*D2, U>, <, , D> 

1400 PLOT <, , D> ! PLOT <, ,U> : GOTO 1430 
1410 GOSUB 1620 ! UJ--W+1 : PLOT <,,U>:GOTO 1430 
1420 PLOT <F1*D1, U3, U>, <, , li>,<, , >, <, , U>! 

E3-INT ( U3) ; E7-INT (F1*D1 ) I GOTO 1570 
1430 E3«INT<F2#D2) ! E7«INT<F1*D1 ) : RETURN 
1440 E1«F2*Y3-INT. < F2*Y3 ) 

1450 E-E1+E4 

1460 P8-F2*D2+E’. P9~INT(P8) 

1 4 70 E5=F 1 *X3- INK F 1 #X3 > 

1480 E6=ES+E8 

1490 S8=F1»D1+E6: S9«INT(S8) 

1500 IF U2~0 THEN 1510: PLOT <, , U>, <S9, 0, U>! GOTO 1530 

1510 IF U3—1 THEN 1520! PLOT < , , D> , <89 , P9 , D> : GOTO 1530 

1520 PLOT <, , U>, <89, P9, U>: U3-0 

1530 IF U2-1 THEN 1540; -E3«€3+P9 

1540 E4«F2*Y3~E3 

1550 E7-S9+E7 

1560 E8=F r l*X3~E7 

1570 RETURN 

1580 DEFFN ' 125: IF ( Y--C2) /Y0+YK9. 999THEN 1590: 02=1 
1590 IF < Y~C2 > /Y0+ Y 1 > .001 THEN 1600; U2*l 
1600 IF < X --01 ) /XO+X 1 <9 » 9 99 THEN 1610! U2-1 
1610 RETURN 

1620 ON (VAL(P$)+1 ) GOSUB 1630, 1640, 1650 ! RETURN 

1630 PLOT <0, 8, U>, <-7, -12, D>, <14,0, D>, <~7, 12, D>, <0, -8, U> i RETURN 
1 640 PLOT <0, 6, U>, <-4,-12, D>, <10,8, D>, <-12, 0, D>, <10, ~8,D>, <-4, 1 
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TABLE VIII CONCLUDED 


D 

16 SO PL 


y • y : ■ 0 » 


•■••6 » U > : 

or <~4 


<0 , -8 , D> 


PI:- TURN 

y 7 y 0 y ***. 8 y 


y < 1 6 y By D > y <.Q y 

1 660 Dlr-FFN ' 06; SELECT 
1670 PRINT "POSITION 


0 y C) > y • •• S y 1 6 y D > , < 8 y 0 , D > « <•• ■ 4 y 7 y D > y < -• 7,4, D > , 
8 y D > , < - 7 , 4 y 0 > y < 0 , 0 , U > I R E r U R N 
:;i RlNT OOS; PRINT HEX (03); PLOT <i,,C>,<,.S> 
LABEL WHERE IT WILL NOT INTERFERE WITH CRAP 


1680 P$* " " : KEY IN P* , 1 700, 1 7 1 0 
1690 GOTO 1680 

1700 IF P$~HEX ( OD ) THEN 1780; IF P$“HEX<02) THEN 1730: IF P$«HEX(0 
8) THEN 1740! PLOT <, , P$>, <13, , U>i GOTO 1680 

1710 IF P$~-HEX < OS ) THEN 17S0: IF P*«HEX<06> THEN 1760; GOTO 1770 

1720 PLOT «), --20 , U>, <--999,0, IJ>: GOTO 1680 

1730 PLOT <13y0,U>; GOTO 1680 

1740 PLOT <™13,0,U>; GOTO 1680 

1750 PLOT <0y20yU>; GOTO 1680 

1760 PLOT <0,-20yU>: GOTO 1680 

1770 STOP 

1780 RETURN 

1790 PLOT <1, , C>, <14, y 8>y <y y G$> 

1800 PLOT < 1 1 2 y - 2 S y U > ! P LOT 
1810 PLOT <--238, -2SyU>; PLOT <,,•!$> 

1820 PLOT <--140, -2£yU>! PLOT < n R$> 

1830 PLOT <--98y--2SyU>: PLOT <„ » G$> 

1840 PLOT <,,R> 

I860 RETURN 


1860 PRINT "LOAD DATA TAPE ON UNIT B INPUT "FILE l! UMBER TO BE 
LOADED" yNl : SELECT TAPE 1 OB : P9™0; REWIND IFOR K=1 TO 50;P9«P9+l! 

NEXT K : IF N 1 « 1 THEN 70; N=N1~1 ; SKIP NF 
1870 DATA LOAD "LTEST" 1 DATA LOAD D$ , , Z1 „ D < 8 ) y Q ( ) , H ( ) r P ( ) , G1 ( ) ; 


IF END THEN 1880; GOTO 1880 
1880 SELECT TAPE 10A ; K6“D ( 8 ) : X$=" TIME IN MINUTES" 
I..AT ION IN LBS" 

1890 RETURN 


Y$“" C02 ACCUMU 


BLOCKF COMPUTER PROGRAM 


Pile Name "BLOCKF" 


Abstract This program will numerically approximate the radiation 

view factor (F) between a rectangle and an N-sided poly- 
gon. Internally, the program subdivides the N-sided 
polygon into N-2 triangles and calculates the view factor 
from the rectangle to the N-2 triangles. With the use 
of form factor algebra, the program combines the triangles 
back into a polygon. This program only operates on two 
surfaces at a time and will not include the effects of a 
shadowing surface. For a detailed explanation of a ra- 
diation view factor, see reference 1. The program is de- 
signed for use with a Wang 2200 series mini computer. A 
sample case is included in figure 10. 


Program Description 


The program numerically calculates 
the following equation: 


F Al “ A 2 


1 

*1 


Jki Ja 2 


the view factor (F) according to 

C0S©1 COS 0 2 dA X dA 2 
7TR 2 


?! ?3 



The required inputs to perform the calculations are as follows: 

1. Number of sides of A 2 (the polygon, MAX=20) 

2. X,Y,Z co-ordinates of each corner of the polygon (must be 
entered sequentially counter clockwise) . 
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3. XYZ of three corners of A]_, the rectangle (must be entered 
sequentially counter clockwise) . 

4. Grid size of the rectangle. The grid size is actually the 
size of the element dAi. The more subdivisions the more 
accurate the results and more computer time used. 

The program outputs the following results on the CRT screen: 

1. The area of each surface, 

2. The view factors Fq_2 and F2-1* 

At this point, the user is given the option to select a different grid 
size (see sample case) . 

A sample case of 'BLOCKF" is provided in figure 11 and a listing of the 
program is provided in Table IX. 


Reference 1 

Hamilton D.C. and W. R. Morgan, "Radiant Interchange Configuration 
Factors," NACA Technical Note 2836, 1952. 
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FIGURE 11 

BLOCKF SAMPLE CASE 



INPUTS 


OUTPUTS 


x y z 

0 0 0 

10 0 

110 


GRID SIZE 
2 
4 
8 

16 


F 

.198 

.19956 

.19993 

.20001 


x y z 

0 0 0 

0 10 

Oil 
0 0 1 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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TABLE IX BLOCKF PROGRAM LISTING 


10 
20 
30 
40 
SO 
SO 
70 
80 
90 ! 
100 
110 
120 
130 
140 
ISO 
ISO 
170 
180 
190 
200 
210 
220 
230 
240 
2 SO 
2 S 0 
270 
280 
290 
300 
310 
320 
330 
340 
3 S 0 
3 S 0 
370 
380 
390 
400 
4 1 0 
420 
430 
440 
4 SO 
460 
470 
480 
490 


D X M X ( 2 1 > „ Y < 2 1 ) „ Z C 2 1 > . U ( 3 > „ V < 3 ) ,, W C 3 ) 


INPUT 11 NUMBER OF SIDES OP VIEWED POLYGON. 

PRINT "INPUT PI PM OF VIEWED POLYGON " 

FOR 1= 1 TO M 

PR I NT " P " 1 I ; : INPUT X(I5,Y(I),Z( I ) 

NEXT I 

X < M+' 1 ) “X <10! Y C M+.L ) =Y ( 1 ) ! Z ( M+ 1 ) «Z < 1 ) 

PRINT " INPUT P1,P2,P3 OF VIEWING RECTANGLE 
FOR 1=1 TO 3 

PR I NT 11 P " ! I ; ! I NPUT LJ ( I ) , V ( I ) „ W < I ) 

NEXT I 

PRINT : PRINT 
FOR 1=1 TO M 
PRINT "P" ; I ; " 

NEXT I 

STOP " KEY CONTINUE " 

PRINT ; PRINT " VIEWING SURFACE 
FOR 1=1 TO 3 

PR I NT " P " I " = " ; u < I ) ; " , " ; V <1 ) ; " f " ; W < I ) 

NEXT I 

M2=0:.FQR 1=2 TO <M-1) 

A1“X ( 1 ) ••••X ( I ) 1 A2=Y ( 1 ) -Y < I ) ! A3»Z < 1 ) -Z ( I ) 

B1=X< 1 )--X< 1 + 1 > ; B2=Y< 1 >-Y< 1 + 1 > 1 B3*Z< 1 )-Z( 1 + 1 > 
C 1 = A2*B3- A3*B2 1 C2=A3*B 1 -A 1 «B3 : C3=A1*B2--A2*B 1 
C0=C1#C1 +C2*C2+C3*C3 ! CO«SGR < CO ) 

M2=M2+C0 

S*M2 


M < MAX- 20 ) " , M 


VIEWED SURFACE (SURFACE 2 ) 

•" ; x < I > j %" ; Y < I ) . j " , " ; z < I ) 

(SURFACE 1 ) 


NEXT I 
A 1 “U (3 ) ■ 
A2=V ( 3 ) • 
A3=W ( 3 ) • 
B 1 = U( I >• 
B2=V ( 1 > • 
B 3 =W ( 1 ) ■ 
REM C = 


M 2 - 
•U ( 2 ) 

■ V ( 2 > 

• W ( 2 ) 

•U ( 2 ) 

•V ( 2 ) 

■ W ( 2 ) 

A CROSS 


B 


“ - GO ; G 1 =G 0~1 
VIEWING SURFACE 


C 1 = A 2 * B 3 ™ A 3 * B 2 
C2=A3*B1 ••••«! *B3 
C 3 = A 1 * B 2 A 2 * B 1 

M 1 = ( C 1 «-C 1 +C2*C2+C3*C3 ) i M 1 =SQR ( M 1 ) 

N1-C1/M1-: N2=C2/M1 : N3=C3/M1 
INPUT "GRID SIZE OF VIEWING SURFACE 
PRINT HEX (03) '.PRINT " GRID SIZE Of 
D 1 = IJ ( 2 ) + S )(- C A :!. +B 1 ) / GO 
D2=V ( 2 ) + S« ( A2 »-B2 ) /GO 
D3'*W (2 ) + . S« ( A3+B3 > /GO 

SELECT R ; S9 :::: 0 1 REM CALCULATE LOCATION ON VIEWING SURF AC 
FOR K>0 TO Gl! FOR J=0 TO 01 
R 1 =D 1 + ( J* A 1 ,+K*.B 1 ) /GO 
R2=D2-v ( J»A2+K-XB2 ) /GO 
R3=D3+ ( J-«- A3+|<«B3 > /GO 


" ; GO 


ORIGINAL PAGE IS 
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TABLE IX CONCLUDED 


SOO FOR TO H 
£10 R4*X(L)-R1 
520 R5“Y ( L) -R2 
£30 R6='Z ( L > -R3 
£40 R7=X ( L •*•:!. ) -X ( L ) 

£50 R8~Y ( L-i- 1 } ~ Y < L ) 

£60 R?=Z<L+1 )-Z(L) 

£70 S1=R7* ( R5*N3~E6*N2 ) ~R8* < R4#N3-R6#N1 > +R9* ( R4*N2~RS*N1 ) 

£80 S2=R4*R4+R5#R5+R6*R6 
£90 S3- 2* < R 4 * R 7 + R £ * R 8 R 6 * R 9 ) 

600 S 4 ~ R 7 * R 7 + R 8 * R 8 + R 9 * R 9 
61 0 S£«4*S2#S4-S3*S3 
620 86=SQR<85) 

6 3 0 S 7 " 8 1 / < ti p I *S6) 

640 8 8 ~ A T N < (2*84+83)/ S 6 ) - ATN ( 83/86 > 

650 88-88*87 
660 89=89+88 

670 NEXT L 1 NEXT J : NEXT K 
680 89=89/ ( GO* GO ) 

690 PRINT ; PRINT " F ( 1 , 2 ) = "JABS (89) 

7 0 0 P R I N T " E ( 2 , 1 ) - " ; A B S ( 8 9 ) * M 1 / M 2 

710 PRINT " A ( 1 ) = " ; Ml 

720 PRINT " A ( 2 ) = M2: PRINT 

730 G1«0: INPUT "ANOTHER GRID SIZE (1*YES, 0=N0>%RliIF G1 =0 THE 
N 740: GOTO 400 r 

740 END 
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SCRIPTF COMPUTER PROGRAM 


File Name 


"SCRIPTF" 


Abstract This program calculates radiation exchange factors from 

a user provided view factor matrix, area vector and 
emitance vector. The program is designed for use with 
a Wang 2200 series mini computer system with the Matrix 
ROM Option. 

Program Description 

J 

This program calculates the radiation exchange factors for use in 
thermal network analysis programs. The 5” is defined as follows: 

(T-' 

J I-*J = qNET I-» J 

6A I (Ti^ - T2 4 ) 


For a gray diffuse emitter, diffuse reflector with constant properties 
the & can be obtained closed form. The method of solution most suit- 
able for the Wang 2200 involves the use of matrix algebra. Derivation 
of the matrix formulas is fully explained in Reference 2. 

The program requires as an input the following values: 

1) Number of surfaces. 

2) The area and emissivity of each surface (any units for the 
area may be used). 

3) Half of the view factor matrix (the program calculates the 
other half through form factor algebra) . 

An example case would be the radiation network inside of a cube one 
foot on each side with all surfaces having an emitance of 0.5. The 
view factor between sides of a cube at 90 degrees to each other was 
calculated in the sample case for the "BLOCKF" program and is .20004. 
The view factor between opposing sides is .19984. Therefore, with 
1 and 2 as the top and bottom and 3 through 6 as the sides of the cube, 
the view factor matrix would be as follows: 

I J VIEW FACTOR 


0 

.19984 

.20004 

.20004 

.20004 

.20004 

0 

.20004 

.20004 

.20004 

0 

.20004 

.19984 

.20004 


VIEW FACTOR 


I J 


4 4 0 

4 5 .20004 

4 6 .19984 

5 5 0 

5 6 .20004 

6 6 0 

Figure 12 is the program output for the sample case. The program first 
prints out an image of the input data and the full view factor matrix. 
Next it prints out the matrix. 

A listing of the "SCRIPTF" program is included in Table X. 

Reference 2 

Wiebelt, John A., Engineering Radiation Heat Transfer , August, 1966, 
pages 115-121. 


FIGURE 12 


SCRIPTF SAMPLE CASE 


SURF # 

AREA 

EMISS 

1 

1 

.. 5 

pi 

1 

A 

u 

5 

1 

5 

4 

1 

. 5 

1 

1 

" § 


F 

MATRIX 


F ( 

1 

,, 1 ) 

0 

F< 

1 

, 5 ) 

„ .19384 

F ( 

1 

, 3 ) 

. 20004 

F ( 

1 

, *'+ > 

. 20004 

F ( 

1 

, E> ? 

u 20004 

F ( 

I 

6 > 

. 20004 

SUM 

1 

i 

F*iOO 

too 


F t 2 , 1 5 

« 19984 

F< 2 , 2 ) 

0 

F C 2 , 3 ) 

. 20004 

F( 2 , 4 ) 

. 20004 

F ( 2 , 5 ) 

. 20004 

F< 2 ‘ 6 ) 

„ 20004 

SUM F*I00 

100 


F ( 

3 , 1 ) 

u 20004 

F ( 

3,2) 

„ 20004 

F < 

3,3) 

0 

F < 

3,4) 

, 20004 

F< 

3,5) 

. 19984 

F< 

3,6) 

. 20004 

SUM 

F*100 

100 


F < 

4,1) 

. 20004 

F < 

4,2) 

. 20004 

F ( 

4,3) 

. 20004 

F ( 

4 , 4 ) 

0 

F ( 

4 , 5 ) 

. 20004 

F 

4 , 6 ) 

. 13984 

SUM F* 1 00 

100 
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F ( 

S 

*1 

1 

) 

. 20004 

P ( 

5 

9 

2 

) 

. 20004 

F ( 

S 

9 

3 

) 

. 1 9984 

FC 

5 

9 

4 

) 

. 20004 

F ( 

5 

3 

s 

j 

0 

F< 

5 

3 

6 

) 

. 20004 

SUM F*1QG 


.too 


F( 

6 , 

1 ) 

. 20004 

F ( 

8 , 

2 ) 

. 20004 

F ( 

6 , 

3 ) 

. 20004 

F< 

8 „ 

4 j 

. 19984 

F < 

6 , 

S ) 

. 20004 

F ( 

8 j 

6 ) 

0 

SUM F*100 

too 


SCRIPT F MATRIX 


SF( 

1 

9 

1 

) 

4. 54545484E-02 

SF( 

1 

9 

2 

) 

9. 08780338E-02 

SF( 

1 

9 

3 

) 

9. 09173544E-02 

SFt 

1 

9 

4 

) 

9. 09173544E-02 

SF ( 

1 

9 

s 

) 

9. 09 t 73544E-02 

SF ( 

1 

P 

8 

) 

9. 09173544E-02 

SUM 

SF/e 

’ TC 

too 

99. 99999999978 


SF ( 

2 , 

1' ) 

9. 08760338E-02 

SF ( 

2 ? 

2 ) 

4. 54545484E-02 

SF( 

c. ? 

3 ) 

9. 09173544E-02 

SFC 

2 , 

4 ) 

9.09173544E-02 

SF < 

• ^ j 

5 ) 

9. 09173S44E-02 

SF ( 

2 , 

6 ) 

9 . 09 173544E-02 

SUM 

SF/E 

I *100 

99. 99999999978 
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SF ( 

3,1) 

SF( 

3,2) 

SF ( 

3,3) 

SF( 

3 „ 4 ) 

SF < 

3,5) 

SF ( 

3 ? 6 ) 

SUM 

SF/E *100 


SF ( 

4 j) 

1 

) 

SF ( 

4 ? 

s 

) 

SF ( 

4 , 


) 

SF C 

4 ? 

4 

) 

SF 

4 , 

5 

) 

SF" ( 

4 

6 

) 

SUM 

SF/E 

- *8 

100 


SF( 

5 

3 

1 ) 

SF C 

S 

5 

2 ) 

SF ( 

5 

3 

3 ) 

SF < 

5 

3 

4 ) 

SF < 

5 

3 

5 > 

SF C 

5 

3 

6 ) 

SUM 

SF 

/ 

/ 

E *100 


SF C 

G 


1 

) 

SF ( 

6 

3 

2 

) 

SF C 

6 

3 

3 

) 

SF ( 

G 

3 

*4* 

) 

SF ( 

G 

3 

5 

) 

SF C 

G 

3 

6 

) 


9. 09173544E-02 
9, 09173544E--02 
4.5454S484E-02 
9. 09173544E-0E 
9. 0876033SE-02 
9.09173544E-G2 
99. 9999999998 


9 . 09 1 73544E-02 
9.09173544E-02 
9.09173544E-02 
4. 54545484E-02 
9. 09173544E-02 
9. 087G0336E-02 
99.99999999978 


9. 09173544E-02 
9. 09173544E-02 
9. 087 603 36E-02 
9 . 09 1 73544E-02 
4. B4545484E-02 
9. 09173544E-02 
99. 9999999998 


9. 09173544E-02 
9. 09173544E-02 
9. 09173544E-02 
9. 0S760336E-02 
9.09173544E-02 
4. S4545484E-02 


SUM SF/E *100 


99. 99999999978 
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TABLE X 


SCRIPTF PROGRAM LISTING 


10 DIM A < 24 > , E ( £4 > , B ( 24 ) , S ( 54 > , F ( 24 , 24 > , C ( 24 , 24 ) : N9=24 
20 INPUT "NUMBER OF SURFACES" , Ns IF N>N9 THEN 20 
BO MAT REDIM A (N) „ E(N) , B ( N ) , S ( N) „ F(N„ N) „ C < N, N } 

40 FOR 1=1 TO N 

50 PR I NT " SURFACE " S I S S INPUT " AREA, EMISS “ , A ( I ) , E(I ) 

GO NEXT I 

70 FOR 1 = 1 TO NS FOR J=I TO N 

SO PRINT "Ft " IS", " ; Js " ) " ; S INPUT FC I , J) 

90 F < J , I)=F(I, J)*A(I)/A(J) 

100 NEXT J S NEXT I 

110 SELECT PRINT 215 ( 120) : PRINT SPRINT 

120 PRINT "SURF # AREA EMISS" 

130 FOR 1=1 TO N: PRINT I , A(I ) , E( I ) S NEXT I: PRINT SPRINT 
140 PRINT " F MATRIX" SPRINT SPRINT 
150 FOR 1 = 1 TO NS E7=0S FOR J=1 TD N 
160 E7=E7+F(I,J> 

170 PRINT "F(": IS ", ";JS" ) " ; F ( I , J > 

ISO NEXT -Js PRINT " SUM F*100 " SE7*100 

190 PRINT SPRINT S NEXT I 

200 PRINT SPRINT SPRINT SPRINT " SCRIPT F MATRIX" S PRINT SPRINT 


210 FOR K=1 TO N 
220 FOR 1=1 TO N 
230 B(.I )=-A( I >*F< I,K)*E(K> 

240 FOR J=1 TO N 
250 D= 1 

2G0 IF I=J THEN 230 
270 D=0 

280 C< I, J)=CF< I, J>*(1~E< J) )-D)*A(K)/A( J)/E( J)*A( I) 
290 NEXT J S NEXT I 
300 MAT C = INV <C) 

310 MAT 3=C*B 
320 E7=0 

330 FOR J=1 TO N 

340 PRINT "SFC'SKS ", "S J; » ) ”;S(J) 

350 E7=E7+S(J) 

360 NEXT J 

370 PRINT "SUM SF/E *100 " ; E7/E(K )*100 

380 PRINT SPRINT 
390 NEXT K 

400 SELECT PRINT 005(64) 

410 END 



APPENDIX 1 


DERIVATION OF CO2 REMOVAL RATE FOR PROGRAM "LIOHCRT" 


The test setup is shown in Figure 13 . The partial pressure of CO2 
is known at the ARS inlet and remains constant up to the LiOH canister 
inlets. Therefore, the problem reduces to the configuration shown in 
Figure 14 . 

The LiOH chemical reaction is as follows: 

2LiOH + CO2 (9) = Li 2 CO3 + H 2 0 (9) 

Therefore, for each mole of C0 2 gas removed, approximately a mole of 
H2O vapor is added to the gas stream. This is an approximation be- 
cause a small amount of water vapor collects in the LiOH bed. There- 
fore, 


ni n2 

The ideal gas law shows: 


V = nMR 

T P 

R = R 

M 

V = nR 

u t 

T P 

or ORIGINAL PAGE IS 

OF POOR 4GiU.UA 

Vl - SIR 

Ti Pi 

V2 = H2^ 

T2 ^2“ 

The total pressure across the LiOH bed remains approximately constant, 
therefore, 



Since the number of moles, the total pressure, and remain constant 
across the LiOH bed, 

V]_ v 2 

Ti T 2 


UPPER LiOH 
CANISTER OUTLET 
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Developing the formula for CO 2 mass removal rate Is as follows: 

AMCO 2 (removed) = MC02 (at 1) - MC02 (at 2) 

= p C0 2 l V j - PC0 2 2 V2 
RCO2 Tx RCO2 T2 

but !i = ^2 

T 1 t 2 

. AMCO 2 = (PCO 21 - PC0 22 ) 

rt 2 


RC0 2 = 35.1 —■ t ft 

lb m - or 


V is in Ft 3 /MIN 


PC0 2 is in mmHg 
Ft3 

AMC0 2 = V( — ) -AP (mmHg) • (2.785 lbf/Ft 3 ) 

mmHg 

35.1 V lbm-°F^ (T+460) (°R) 

AMC0 2 = 0.07931 * V * A P/ (T+460) in lbm/min. 


ORIGINAL 
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